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Abstract 

We  are  concerned  here  with  the  nonlinear  one -dimensional  hydromagnetic  flow 
problem  of  the  resolution  of  an  initial  shear  flow  discontinuity  in  a  perfectly 
conducting,  electrically  neutral,  compressible  fluid.  By  one-dimensional  hydro - 
magnetic  flow  we  mean  a  hydromagnetic  flow  in  which  the  dependent  variables,  the 
magnetic  intensity,  fluid  velocity,  density,  pressure  and  entropy  per  unit  mass, 
are  functions  of  one  space  variable,  x  say,  and  time  t;  no  restriction  need  be 
imposed  as  to  the  presence  or  absence  of  the  y-  and  z-components  of  the  magnetic 
intensity  and  fluid  velocity.  Consider  a  plane  surface  normal  to  the  x-axis 
across  which  a  component,  the  y-component  say,  of  the  velocity  is  discontinuous 
and  across  which  no  fluid  passes.  In  addition,  let  it  be  assumed  that  these 
variables  are  independent  of  x  initially  and,  moreover,  that  the  transverse 
component  of  the  magnetic  field  (i.e.,  the  component  normal  to  the  x-axis), 
vanishes  initially.  Such  a  flow  discontinuity  is  referred  to  as  a  shear  flow 
discontinuity.  Since  the  basic  hydromagnetic  equations  are  invariant  under 
galilean  transformations  we  may,  without  loss  of  generality,  assume  that  the 
longitudinal  component  of  the  velocity  (i.e.,  the  component  parallel  to  the  x- 
axis),  and  the  z-component  of  the  velocity  also  vanish  and  that  the  y-component 
of  the  velocity  in  the  region  x  >  0  is  equal  and  opposite  to  that  in  the  region 
x  <  0»   Hereafter,  we  shall  assume  this  to  be  the  case.  Now  if  the  longitudinal 
component  of  the  magnetic  field,  H  ,  is  assumed  not  to  vanish,  it  can  be  shown 
that  a  shear  flow  discontinuity  is  not  an  admissible  hydromagnetic  contact  dis- 
continuity. The  existence  of  an  initial  shear  flow  discontinuity  in  the  presence 
of  a  longitudinal  magnetic  field  must  therefore  be  accommodated  by  a  suitable 
motion  of  the  medium.  We  refer  to  this  motion  as  being  a  'resolution'  of  the 
initial  shear  flow  discontinuity.  The  object  of  this  worK  is  to  construct  the 
solution  of  this  nonlinear  initial- value  problem,  to  discuss  some  of  the  general 
features  of  this  solution  and  to  provide  for  each  time  t  >  0  the  numerical  basis 
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for  obtaining  the  waveforms  of  the  variables  which  characterize  the  flow. 

An  idealized  physical  situation  in  which  a  shear  flow  discontinuity  occurs 
is  the  following.     Consider  a  layer  of  perfectly  conducting,  compressible  fluid 
with  faces  transverse  to  the  x-direction.     Assume,  at  the  initial  instant,   that 
the   fluid  within  the  layer  moves  uniformly  in  the  negative  y-direction,  that  only 
the  fluid  within  the  layer  is  in  motion,  and  that  the  pressure  and  density  are 
continuous  across  the  layer  faces  and  ars  independent  of  x.     Such  an  initial 
motion  may  be  regarded  as  an  idealization  of  a  jet  proceeding  in  the  negative 
y-direction  and  much  wider  in  the  z-direction  than  in  the  x-direction.     At  each 
face  of  the  jet  we  encounter  a  shear  flow  discontinuity.     The  waves  resulting 
from  the  resolution  of  the  shear  flow  discontinuity  at  each  face  will  interact 
only  after  a  certain  time.     Before  this  time  the  situation  may  be  described  by 
a  wave  motion  resulting  from  a  single  interface.     It  is  this  motion  that  is  con- 
sidered here. 

The  methods  employed  to  determine  the  motion  parallel  those  employed  in  one- 
dimensional  gas  dynamical  problems.     It  is  shown  that  the  initial  discontinuity 
is  resolved  by  two  fast  shocks  moving  away  from  the  position  of  the  initial  dis- 
continuity followed  by  slow,   centered  hydromagnetic  rarefaction  waves.     These 
rarefaction  waves  connect  the  constant  state  immediately  behind  the  moving  shock 
fronts  either  to  the  endpoints  of  an  expanding  zone  in  which  the   thermodynamic 
state  is  constant  and  the  medium  is  motionless  or  to  the  endpoints  of  a  vacuum 
zone  which  expands  with  an  appropriate   'escape  speed'.     The  fast  shocks  in  question 
are  fast  gas  dynamical  shocks  unless  the  square  ratio  between  the  sound  the  Alfven 
disturbance  speed  in  the  undisturbed  medium  is  less  than  unity  and  the  shock 
strength  is  less  than  a  certain   'critical'   value.     In  this  last  case  the  resolution 
is  effected  by  switch-on  shocks.     Fast  gas  dynamical  shocks  are  gas  shocks  in  which 
the  speed  of  the  shock  front  is  greater  than  the  Alfven  characteristic  speed  behind 
the  6hock.     Switch-on  shocks  are   shocks  having  a  vanishing  transverse  magnetic  field 
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ahead  of  the  shock  front  but  nevertheless  a  non-vanishing  transverse  ragnetic 
field  component  behind  the  shock. 

Switch-on  shocks,  fast  gas  shocks,  and  slow  simple  wave  solutions,  the 
basic  building  blocks  entering  into  the  construction  of  the  solution,  are  in- 
vestigated in  detail.  It  is  then  shown  that  a  complete  quantitative  description 
of  the  motion  may  be  obtained  from  the  knowledge  of  how  two  pressure  ratios,  p,/p 
and  p,/p  ,  depend  on  the  magnitude  of  the  initial  discontinuity  2u  Q,  the  sound 
speed,  and  the  Aifven  disturbance  speed,  in  the  medium  in  front  of  the  advancing 
shocks.  Here,  p  and  p,  are  the  pressures  in  front  of  and  behind  the  shocks  and 
p_  is  the  pressure  at  the  rails  of  the   simple  waves. 

A  result  of  particular  interest  is  that  in  the  region  between  the  rarefaction 
waves  a  transverse  component,  H     ,,  of  the  magnetic  field  is  created,  even  though 
such  a  component  was  nowhere  present  initially.     When  the  magnitude  of  the  initial 
discontinuity  of  the  velocity  2u     _  is  small,  it  is  shown  that  H     ,  is  proportional 
to  the  product  of  u        and  the   square  root  of  p  ,  the  density  in  the  undisturbed 
medium.     When  u     -  is  large  -  that  is,  when  the  nonlinear  terms  in  the  basic 
equations  are  operative,  it  is  shown  that  H     ,  is  proportional  to   the  product  of 
the   square  root  of  u     0,of  the  fixed  longitudinal  magnetic  field  H  ,   and  of 
the  fourth  root  of  tne  density  p     in  front  of  the  waves.     Thus,  even  when  H     is 
very  small,   a  large  transverse  magnetic  field  is  created  behind  the  rarefaction 
waves  if  the  magnitude  of  the  initial  shear  flow  discontinuity  2u         is  large 

y»u 

enough.     Associated  with  this  result  is  the  fact  that  the  motion  of  the  medium 

differs  radically  from  the  corresponding  gas  dynamical  steady-state   shear  flow 

motion,  even  when  H     is   small,  if  2u     _  is  large  enough.     As  already  indicated 

above,  a  longitudinal  motion  results  while  no  such  motion  occurs  in  the  gas 

dynamical  case.     Moreover,  this  longitudinal  motion  is  such  that  cavitation  is 

produced  whenever  u     „  exceeds  a  certain  minimum  value  which  depends  on  the  ratio 
*  y,0 

of  the  sound  speed  and  the  Aifven  disturbance  speed  in  the  medium  in  front  of  the 
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advancing  shocks. 

The  exact  nature  of  the  dependence  of  p,/p.,  P,/Pn>   H     ../H     and  other  quanti- 

J       X  X       O  J  9  J       X 

ties  of  interest  on  the  parameters  which  characterize  the   state  in  front  of  the 
advancing  shocks  is  obtained  by  a  combination  of  graphical  and  numerical  methods. 
In  addition,   approximate  expressions  for  these  quantities  are  developed  for  small 
and  large  values  of  the  initial  discontinuity.     From  an  analysis  of  these  results 
we  then  ascertain  some  of  the   general  features  of  the  motion,, 
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1.       Introduction 

We  are  concerned  here  with  the  nonlinear  one-dimensional  hydromagnetic  flow 
problem  of  the  resolution  of  a  shear  flow  discontinuity  in  a  perfectly  conducting, 
electrically  neutral,   compressible  fluid.     By  one-dimensional  hydromagnetic  flow 
we  mean  a  hydromagnetic  flow  in  which  the  dependent  variables,   the  magnetic  in- 
tensity, fluid  velocity,  density,  pressure  and  entropy  per  unit  mass,   are  functions 
of  one   space  variable,  x  say,  and  the  time  tj  no  restriction  need  be  imposed  as  to 
the  presence  or  absence  of  the  x  and  z-components  of  the  magnetic  intensity  and 
fluid  velocity.  Because   the  flow  depends  only  on  the  single   space  variable  x,  it  is 
clear  that  shock  and  contact  discontinuities,  when  present,   appear  only  across  planes 
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FIG.  1.  Symmetrical  shear  flow  discontinuity  in  the  transverse  velocity, 

u  .  L  is  the  planar  surface  x  =  0,  which  separates  the  regions  of  equal 

tr 
and  opposite  transverse  flow.  The  quantity  pQ  is  the  density  of  the  fluid, 

p  the  pressure  of  the  fluid,  and  H  >  0  the  fixed  longitudinal  component 

of  the  magnetic  field.  These  quantities  are  assumed  to  be  continuous  across 

L.  The  transverse  magnetic  field  H,   and  the  longitudinal  fluid  velocity  is 

assumed  to  be  zero  throughout.  The  magnitude  of  the  initial  discontinuity 

is  2u  n.  Unsymmetrical  shear  flow  discontinuities  in  u.   may  be  reduced 

to  the  present  symmetrical  case  by  means  of  a  galilean  transformation. 
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normal  to  the  x-axis.  Consider  a  plane  surface  across  which  a  component,  the  y- 
component  say,  of  the  transverse  velocity  is  discontinuous  and  across  which  no 
fluid  passes.  In  addition,  let  it  be  assumed  that  the  remaining  variables  are 
continuous  across  the  surface,  that  these  variables  are  independent  of  x  initially, 
and  moreover,  that  the  transverse  component  of  the  magnetic  field  as  well  as  the 
longitudinal  fluid  velocity  vanish  initially  -  hare  the  longitudinal  and  trans- 
verse components  of  any  vector  quantity  are  the  components  of  this  quantity  parallel 
and  perpendicular  to  the  x-axis,  respectively.  Such  a  flow  discontinuity  will  be 
referred  to  as  a  shear  flow  discontinuity.  If  the  component  of  the  magnetic  field 
normal  to  the  discontinuity  surface  vanishes, the  discontinuity  reduces  to  a  con- 
tact discontinuity  in  one -dimensional  gas  dynamical  flow.  As  is  well  known,  such 
a  discontinuity  may  persist,  theoretically,  for  all  time.  If,  on  the  other  hand, 
the  longitudinal  component  of  the  magnetic  field  does  not  vanish,  it  can  be  shown 
that  this  discontinuity  is  not  an  admissible  hydromagnetic  contact  discontinuity* 
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FIG.  2.     Jet  proceeding  in  the  negative  y-direction.     At  each  face  of  the 
jet  we  have  a  shear  flow  discontinuity. 
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The  existence  of  an  initial  shear  flow  discontinuity  in  the  presence  of  a  longi- 
tudinal magnetic  field  must  therefore  be  accommodated  by  a  suitable  motion  of  the 
medium.  We  refer  to  this  motion  as  being  a  'resolution'  of  the  initial  shear 
flow  discontinuity,  or  as  the  motion  which  'resolves'  the  initial  shear  flow  dis- 
continuity. The  object  of  the  present  work  is  to  construct  such  a  resolution,  to 
describe  the  general  features  of  this  motion,  and  to  provide,  for  each  t  >  0,  the 
numerical  basis  for  obtaining  the  waveforms  of  the  variables  which  characterize 
the  flow.  Here  and  in  subsequent  sections  we  may  assume,  without  loss  or  gen- 
erality, that  the  transverse  velocity  on  one  side  of  the  discontinuity  surface 
is  equal  and  opposite  to  that  on  the  other  side  and  that  it  is  parallel  to  the 
y-axis.  In  Figure  1  we  have  depicted  one  such  symmetrical  initial  situation. 
Non-symmetrical  shear  flow  discontinuities  may  be  reduced  to  the  symmetrical 
cas  by  means  of  a  galilean  transformation.  In  Figure  3  we  have  sketched  the 
waveforms  of  the  motion  corresponding  to  a  shear  flow  discontinuity  of  the  type 
shown  in  Figure  1. 

An  idealized  physical  situation  in  which  a  shear  flow  discontinuity  occurs 
is  the  following.  Consider  a  layer  of  perfectly  conducting,  compressible  fluid 
with  faces  transverse  to  the  x-direction  [see  Figure  2j.  Assume,  at  the  initial 
instant,  that  the  fluid  within  the  layer  moves  uniformly  in  the  negative  y-direction, 
that  only  the  fluid  in  the  layer  is  in  motion,  and  that  the  pressure  and  density 
are  continuous  across  the  layer  faces  and  are  independent  of  x.  Such  an  initial 
motion  may  be  regarded  as  an  idealization  of  a  jet  proceeding  in  the  negative  y- 
direction  and  much  wider  in  the  z-direction  than  in  the  x-direction.  A  similar 
situation  was  described  by  AlfWnL  ■•,  in  explaining  the  possible  origin  of  hydro- 
magnetic  wave 8.  At  each  face  of  the  slab  we  encounter  a  shear  flow  discontinuity. 
The  waves  resulting  from  the  resolution  of  the  discontinuities  at  each  face  will 
interact  only  after  a  certain  time.  Before  this  time  the  situation  may  be  des- 
cribed by  a  wave  motion  resulting  from  a  single  interface.  It  is  this  motion  that 

is  considered  here. 

The  methods  employed  to  determine  the  motion  parallel  those  used  in  one-dimen- 
sional gas  dynamical  problems.  In  such  problems  initial  discontinuities  of  a 


.  In- 
sufficiently simple  character  are  resolved  by  Beans  of  shocks  and  simple  waves. 
As  will  be  recalled,  the  usefulness  of  simple  waves  is  a  consequence,  primarily, 
of  the  fact  that  a  flow  region  adjacent  to  a  region  of  constant  state  is  always 
a  simple  wave.     The  corresponding  simple-wave  and  shock-wave  building  blocks  are 
also  available  in  the  theory  of  one -dimensional  hydromagnetic  flow.     The  situation 
here,  however,  is  complicated  by  the  fact  that  there  are  three  possible  types  of 
shocks,  the  so-called  transverse  shocks  and  the  slow  and  fast  shocks,  and  in 
addition,  the  corresponding  types  of  simple  waves.     Nevertheless,  from  an  examina- 
tion of  the  known  properties  of  the  various  shocks  and  simple  waves  one  is  led,  in 
the  present  problem, to  a  resolution  consisting  of  a  combination  of  fast  shocks  and 
slow,  centered,  simple  waves.     Because  of  the  initial  conditions  the  transverse 
components  of  the  magnetic  field  must  vanish  ahead  of  the  shocks.     As  will  be 
seen  later  this  implies  that  the  fast  shocks  must  reduce  either  to  fast  gas 
dynamical  shocks  or  the  so-called  switch-on  shocks.     Switch-on  shocks  are  fast 
shocks  having  a  vanishing  transverse  magnetic  field  ahead  of  the  shock  front  but 
nevertheless  a  non-vanishing  transverse  magnetic  field  behind  the  shock  front. 
Fast  gas  dynamical  shocks  are  gas  shocks  in  which  the  speed  of  the   shock  front  is 
greater  than  the  so-called  Alfve"n  characteristic  speed  behind  the  shock.     In  terms 
of  these  building  blocks  the  resolution  of  the  shear  flow  discontinuity  may  be 
briefly  described  as  follows:     The  fluid  velocity  is  an  odd  function  of  x,   and  the 
magnetic  field,  the  pressure  and  the  density  all  are  even  functions  of  x.     The 
initial  discontinuity  is  resolved  by  two  fast  shocks  moving  away  from  the  position 
of  the  initial  discontinuity  followed  by  two  slow,  centered,  simple  rarefaction  waves. 
These  rarefaction  waves  connect  the  constant  state  immediately  behind  the  moving 
shock  fronts  either  to  the  endpoints  of  an  expanding  zone  in  which  the  thermodynamic 
state  is  constant  and  the  medium  is  motionless,  or  to  the  endpoints  of  a  vacuum  zone 
which  expands  with  an  appropriate  longitudinal  'escape  speed' •     The  latter  alternative 
occurs  whenever  the  initial  discontinuity  in  the  transverse  velocity  exceeds  a  cer- 
tain minimum  value,  the    'cavitation  value*  say,  which  depends  upon  the   state  of  the 
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undisturbed  medium.  In  this  case  the  fluid  velocity  has  in  addition  to  the 
longitudinal  escape  speed  a  non-vanishing  transverse  component  at  the  'tails •  of 
the  rarefaction  waves. 

It  will  be  shown  that  a  complete  quantitative  description  of  the  motion  may 
be  obtained  from  the  knowledge  of  how  two  pressure  ratios  P^P.  and  P^/P-i  depend 
on  the  magnitude  of  the  initial  discontinuity  and  on  the  remaining  state  variables 
in  the  undisturbed  medium.  Here  p  and  p,  are  the  pressures  in  front  of  and  behind 
the  shocks  and  p,  is  the  pressure  at  the  'tails' of  the  simple  waves.  The  dependence 
°f  P-i/P  and  Po/Pi  on  ^ne  parameters  in  question  is  established  by  employing  the 
condition  that  the  fluid  velocity  vanishes  in  the  expanding  zone  between  the  rare- 
faction waves  whenever  the  value  of  the  initial  discontinuity  in  the  transverse 
velocity  is  less  than  the  cavitation  value  mentioned  above,  or  the  condition  that 
the  electric  field  vanishes  at  the  tails  of  the  rarefaction  waves  whenever  the 
initial  discontinuity  in  the  transverse  velocity  is  greater  than  the  cavitation 
value.  The  first  condition  is  simply  a  consequence  of  the  fact  that  the  fluid 
velocity  is  an  odd  function  of  xj  the  second  condition  will  be  shown  to  follow 
from  the  quasistatic  approximation  to  Maxwell's  equations  in  the  vacuum  zone  and 
the  hydromagnetic  analog  of  the  Rankine-Hugoniot  conditions  at  the  boundary  of  a 
perfectly  conducting  compressible  medium  which  expands  into  a  vacuum.  In  general, 
it  might  be  expected  that  these  pressure  ratios  depend  on  u  _,  p ,  p  and  H  , 
where  2u  .  is  the  magnitude  of  the  initial  discontinuity  and  p  ,  p  ,  H  are,  in 

y,U  u    u    a 

the  order  given,  the  density,  pressure  and  longitudinal  magnetic  field  in  vhe  un- 
disturbed medium.  We  shall  show,  however,  that  P^/Pq  and  P-$/V±  depend  only  on 

?  ?  2-1  1/2 

the  ratios  s  -  a  /b  and  u  Q/&o   or  alternatively  u  0/bQ.  Here  bo  -  (mHxpq  ) 

is  the  Alfven  speed*  in  the  undisturbed  medium  and  bq   is  the  square  ratio  of  the 

sound  speed  a  »  (yp  p  )  '  and  the  Alfven  speed  in  the  undisturbed  medium.  The 


o"o 


*  The  quantity  \i  is  the  magnetic  specific  inductive  capacity  of  free  space.  The 
Giorgi  mks  system  of  units  is  employed  throughout  this  report. 
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constant  y  is  the  exponent  of  the  density  in  the  adiabatic  equation  of  state  of 

the  medium.  Plots  of  p,/p  versus  u  ./a  and  p-,/p,  versus  u  n/b     are  given 

*r*o        y,0  o    r3  1       y,0'  o    B 

in  Figures  15  and  16. 

As  already  indicated,  a  complete  description  of  the  solution  in  the  various 

flow  regions  may  be  obtained  from  this  basic  information.  Of  particular  interest 

is  the  fact  that  in  the  region  behind  the  rarefaction  waves  a  transverse  component 

of  the  magnetic  field,  H_  ,,  is  created  even  though  such  a  component  was  nowhere 

y>  j 

present  initially.     In  particular  it  can  be  shown  that 


(i.D  «V,3  ~  (p»i/2  Vo 

[cf.    (5.50),   (5.65),  or   (6.10),    (6.1C)]  when  u         is  sufficiently  small  and  that 

i-        -i  -i 

vl/2 


(1.2)        ^3  ~  [{tar)  p./,]1"1  (Vyj0,J 


[cf.  (5.37),  or  (6.11),  (6.11')3  when  u  _.  is  sufficiently  large.  As  already 

y>u 

mentioned,  p  is  the  density  in  the  undisturbed  medium,  H  the  longitudinal 
magnetic  field  ,  \i  the  magnetic  specific  inductive  capacity  of  free  space,  y 
the  adiabatic  exponent,  and  u  0  half  the  initial  discontinuity  in  the  transverse 
velocity.  Thus,  even  when  H  is  very  small,  a  large  transverse  magnetic  field 
is  created  in  the  region  between  the  rarefaction  waves  when  u  n  is  large  enough. 

y>o 

Associated  with  this  result  is  the  fact  that  the  motion  of  the  medium  differs 
radically  from  the  steady  state  gas  dynamical  shear  flow  even  when  H  is  small, 

In  all  of  the  graphs  appearing  in  this  report  the  value  of  the  adiabatic  expon- 
ent y  is  assumed  to  be  5/3.  The  equation  of  state  of  the  medium  is  therefore 
the  same  as  that  of  an  ideal  monotonic  gas. 

In  one-dimensional  hydromagnetic  flow,  the  longitudinal  component  of  the  magnetic 
field  H^  can  be  shown  to  be  independent  of  both  the  space  variable  x,  and  the  time 
variable  t. 
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provided  u    „  is  large  enough.     In  fact  cavitation  will  result  when  u    rt  is 

r  y,0  7,0 

greater  than  the  cavitation  value  noted  above.     The  detailed  nature  of  the  de- 
pendence of  H     ,  on  u    0  is  given  in    Figure     17a     where  we  have  plotted  H     ~/H 

versus  u     „/a     and  in     Figure     17b    where  we  have  plotted  H     0/H     versus  u    ,-,/b^. 
y,0     o  y,3     x  y,u     o 

In  addition,  we  have  plotted   (u     ,)       /b       u     _/b     and  p,A>     versus  u    rt/b    ^fsee 
■  x,3  esc     O)    y,3     o  r3     o  y»0     o   •- 

Figures  18  and  191.     Here,   (u     ,)         is  the  longitudinal  escape  speed  and  u     , 

x,.?  esc  y,j 

the  transverse  fluid  speed  at  the  tail  of  the  rarefaction  wave  in  the  region  x  >  0. 
The  interpretation  of  these  graphical  results  is  deferred  to  Section  6. 

Hydromagnetic  shocks  have  been  treated  by  De  Hoffman  and  Teller  (1950)  Ll, 
by  Heifer  (1952)^   and  Lust  in  (1953)  M"(a)  and  ^  (1955)  M"(b)  #     A  detailed 
treatment  of  gas  dynamical  shocks  may  be  found  in  the  text  of  Courant  and 
Friedrichs'-?-'.     a  theory  of  shocks  and  simple  waves  for  general  systems  of  hyper- 
bolic partial  differential  equations  of  the  type  considered  in  this  report  has 
been  given  by  Lax  (195U)  ^  •*•     The  problem  of  the  resolution  of  a  shear  flow  dis- 
continuity was  first  formulated  by  Friedrichs  (195U)^-     .     This  formulation  was 
based  on  a  new  survey  of  hydromagnetic  shocks  in  which,   among  other  results,  the 
inadmissibility  of  a  shear  flow  discontinuity  in  the  presence  of  a  longitudinal 
magnetic  field  was  noted  and  in  which  some  properties  of  switch-on  shocks  were 
given.     Friedrichs  also  discussed  simple  wave   solutions  of  the  one-dimensional 
system  of  hydromagnetic  equations  treated  here.     Of  particular  interest  is  the 
fact  that  he  succeeded  in  reducing  the  problem  of  obtaining  explicit  slow  and 
fast  simple  wave   solutions  in  media  governed  by  the  adiabatic  equation  of  state 
to  essentially  that  of  solving  a  single  linear  ordinary  differential  equation. 
The  present  work  takes  as  a  point  of  departure  the  formulation  and  results  of 
Friedrichs  just  noted.     In  addition  to  giving  a  detailed  treatment  of  the   shear 
flow  problem  we  extend  slightly  that  part  of  Friedrich's  results  on  switch-on 
shocks  and  simple  waves  which  serves  as  a  basis  for  this  problem. 

The  shear  flow  problem  treated  here  is  perhaps  the  simplest  hydromagnetic 
problem  with  discontinuous  initial  values  that  can  be  formulated.     This  simpli- 
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city  is  clearly  bought  at  the  expense  of  physical  content.     The  loss  of  physical 
content,   however,   is  partly  compensated  by  the  fact  that  a  rigorous  solution  of 
a  nonlinear  problem  is  obtained.     In  addition,  one  may  hope  that,  as  in  gas 
dynamics,   a  familiarity  with  the  properties  of  this  and  similar  one-dimensional 
solutions  may  lead  to  an  increased  understanding  of  more  significant  hydromagnetic 
phenomena. 

In  Section  2  we  give  a  brief  survey  of  pertinent  facts  about  one-dimensional 
hydromagnetic  flow  and  then  formulate   the  problem  in  detail.     Sections  3  and  h  are 
devoted  to  a  discussion  of  the  basic  building  blocks  from  which  the  solution  is 
constructed)  in  Section  3  we  discuss  switch-on  shocks  and  fast  gas  shocks,  and  in 
Section  U,   slow  simple  waves.     These  basic  motions  are  combined  in  Section  5  to 
form  the  solution.     In  this  section  we  also  sketch  our  procedure  for  obtaining 
the   graphical  results  noted  above.     A  survey,  based  on  the  results  of  Section  5> 
of  some  qualitative  features  of  the  motion  is  given  in  Section  6. 

The  reader  who  desires  to  obtain  an  overall  view  of  our  procedure  and  results 
without  entering  into  all  the  detailed  calculations  is  advised  to  read  Section  2, 
subsections  3.1»  U.l»   5.1»   and  Section  6» 
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2.       Survey  of  known  results*     Formulation  of  the  problem 

2.1     The  equations  of  one -dimensional  hydromagne tic  motion 

In  order  to  introduce  notations  and  terminology  employed  throughout 
this  report,   to  motivate  properly  the  formulation  of  the  problem  and  the  choice 
of  building  blocks  entering  into  the   solution,  and  to  prepare  for  the  consider- 
ations of  subsequent  sections,  we   shall  first  summarize  some  pertinent  information 
concerning  one-dimensional  hydrcmagnetic  motion.     Explicitly,  we  shall  give   a 
brief  survey  of  the  possible  types  of  hydromagnetic  shocks  and  simple  waves  and 
recall  some  properties  of  these  motions.     We  shall  then  discuss  in  detail  the 
formulation  of  the  problem.     The   following  exposition  is  based  on  the  work  of 
Friedrichs  ^J  . 

We  begin  with  a  system  of  partial  differential  equations  governing  the 
motion  of  a  perfectly  conducting,  neutral,  non-viscous,  compressible  fluid.     This 
system,  which  is  essentially  a  one-dimensional  version  of  that  employed  by 
Lund  qui  st  ^-  -I ,  is: 

A  H       ■     const   (independent  of  x  and  t) 

3fftr  9fftr      _    ^tr      -      ^x  n 

*1  "3T-  +  ux  "ax"  *  Hx  "ST  +  Htr  aT    "    ° 

^x  ^x 


V        P  3T  +  pux  1ST  +  k  y  *  ^tr>  '2 
3V  *tr  S"tr 


x   ax 


(2.1) 


A3  £    *    k^    '    ° 


\      h  $  *»2*<»  +  7^}  k  (y1^2*  * *  P  *  **>] 


-^|.(ff.u)     -     0 


as  as 

A5  p  aT  +  p  ux  air   "    ° 

A6  e     •     e(ptS),  de     -     TdS  ♦  pdCp"1). 
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All  quantities  appearing  above  are  assumed  to  be  functions  of  the  time,  t,   and 
the  single   space  variable,  x.     H  is  the  magnetic  intensity,  u  the  velocity  of 
the  fluid  particle;   H     and  u     are  the  components  of  these  quantities  in  the  posi- 
tive x-direction  while  H       and  u       are  the  components  of  these  quantities  in  the 
plane  transverse  to  the  x-axis.     The  quantity  ^i  is  the  magnetic  inductive  capacity, 
p  the  density,  p  the  pressure,   T  the  temperature,   e  the  specific  internal  energy 
(that  is,   internal  energy  per  unit  mass),  and  S  the  specific  entropy.     Equation 
A.   is  simply  a  statement  of  the  first  law  of  thermodynamics  for  a  fluid  particle o 
A,-  establishes  the  adiabatic  character  of  the  flow.     Equations  A^  and  Ar  are 
statements  of  the  conservation  of  mass  and  energy,  respectively.     k0       and  A- 
assert  the  conservation  of  momentum  in  the  longitudinal  and  transverse  directions. 
They  are  the  component  equations  of 

p  |£    +  p(u   .   V)u  +  Vp  ♦  uifcc(V  x  i?)     -     0. 

The  term  uflx(V  x  ff)  in  this  equation  is  the  force  exerted  on  the  current  J  -  V  x  H 
by  the  magnetic  field  nH.     Equations  A     and  A,    are  consequences  of  the  relations 

V   •  H     «     0 

||    ♦  V  x  (flxu)     «     0     , 

the  first  of  which  asserts  the  irrotational  character  of  the  magnetic  field.     The 
second  equation  relates  the  time  rate  of  change  of  the  magnetic  induction  off  and 
the  curl  of  the  electric  field 

(2.2)  ST    -     uiT  x  u 

in  the  perfectly  conducting  medium.  The  fact  that  H  is  independent  of  x  and  t 
is  a  special  feature  of  one-dimensional  hydromagnetic  motion. 

It  is  assumed,  here  and  in  what  follows, that  the  fluid  obeys  the  equation  of 
state  for  a  poly  tropic  ideal  pas.    It  then  follows  from  the  adiabatic  character 
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of  the  flow  A-  and  the  first  law  of  thermodynamics  A,  that 


(2.3) 


P  -  K  PY 


e  «  p  p"  /  (y-1) 


where  y  is  the  adiebatic  exponent  of  the  medium  and  K  is  an  arbitrary  function 
of  the  entropy  S. 


2.2  Hyaromagnetic  shocks 

In  a  region  in  which  the  flow  is  continuous,  equations  A^,  A,,  Ar,  and 

A,  imply  A  .     In  this  case  equation  A.    may  be  discarded.     However,  when  a  region 

of  flow  is  divided  by  a  surface  of  discontinuity,  L  -  L(t),  A.  must  be  retained 

and  together  A     -  A,    lead  to  the  hydromagnetic  analog  of  the  RankLne-Kugoniot 
o         h 

conditions  across  L: 

B  [H  ]     -     0 

o  L  nJ 

B1  (mfr    -     Hnu         =     0 

B.  Uu+(p  +  i-ufl")n+  uJInH*        -     0 

(2.U) 

B-3  m    Cf]  "    Q1  1     "     °>  or  equivalently    [m]     -     0 

Bh  t(j-  u2  +  e  ♦  \  \&2V  )  +  un(p  +  \  nH2;  -  n«n(H  •  u)l    -  0. 


Equation  A-  must  be  replaced  by  the  condition,  m[S]  >  0.  This  implies  that 
the  entropy  'behind'  a  shock  is  not  less  than  that  'in  front'  of  the  shock 
( see  next  page ) . 
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Here  H   ,  u     are  FT  •  n  and  u   •  n  respectively  where  the  vector  n  denotes  the 
n       n 

normal  vector  at  any  point  of  the  discontinuity  surface,  L,   and 

(2.?)  [Q]     "     Q±  "  Q0» 

the  jump  in  any  quantity  Q  across  L,  Q.  being  the  value  of  Q  on  the   side  into 
which  n  points.     In  addition  to  have   set 

(2.6)  f  -     p"1 
and 

(2.7)  m     -     p(un  -  Un) 

where  u     and  D     are  the  fluid  and  shock  velocities  in  the  normal  direction.     Note 
n  n 

that  the   flux  through  the   surface,  m,  is  continuous  across  L    |cf.   B^  above],     De- 
pending upon  whether  n  i  0  or  m  •  0  we  speak  of  a  shock  or  of  a  contact  discon- 
tinuity.    Since  we  are  dealing  exclusively  with  one-dimensional  hy dromagne tic 
flow  depending  on  the  space  variable  x  it  is  clear  that  n  is  directed  along  either 
the  positive  or  negative  x-axis.     In  what  follows,  when  the  discontinuity  is  a 
shock  we   shall  employ  the  convention  that  n  points  in  the  direction  in  which  the 
fluid  crosses  the   surface.     Under  these  circumstances  m  is  positive   .     The  region 
into  which  n  points  will  be  referred  to  as  the  region  behind  the  shockj  the  re- 
maining region  will  be  called  the  region  ahead  of  the   shock.     Fluid  therefore 
crosses  the  shcck  wsve  from  the  region  ahead  of  the   shock  to  the  region  behind  the 
shock.     A  shock  is  said  to  be  forward  facing  or  backward  facing  according  as  n  is 
directed  along  the  negative  or  positive  x-directions.     For  a  forward-facing  shock 
the  region  ahead  of  the   shock  lies  to  the  right  of  the  region  behind  the  shockj 
whereas  for  a  backward-facing  shock  the  region  ahead  of  the   shock  lies  to  the 
left  of  the  region  behind  the   shock. 


When  m  -  0  we   assume  that  n"  is  directed  along  the  positive  x-axis. 
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Employing  the   identity 

(2.8)  [PQ]     -     P    [Q]     +      [P]   Q, 
where 

(2.9)  Q    -     £(Q0  ♦  Qx) 

is  the  mean  value  of  any  quantity  having  the  value  Q    behind  and  Q     ahead  of  the 
shock,  we  may  write  equations  B,   -  B_  as  follows: 

C±  mt    QH]     +     l[un]  -  Hn[u]     -     0 

(2.10)  C2  m    £u  ]  +    [p]   n  +   ufT  ff.QFT]   -  ^nL^  1  '  ° 

c3  m  [r]  -   [un]    -    0. 

These  equations  in  conjunction  with  the  similarly  expanded  form  of  B.    may  be  re- 
garded as  a  system  of  seven  linear  homogeneous  equations  in  the  seven  unknowns 
[p],    £f],    [u~]  and  the  jumps  in  the  transverse  components  of  H,     Whenever  non- 
trivial  solutions  of  this  system  exist  the  determinant,   det(C),  must  vanish. 
It  can  be   shown  that 


det(C)     -     m(m2;t    -  txH2.) 


(m2r    -  tiH^Jra2  at  -   (l+a)p  -   Jk)2  +  J^J 


(2«H)  2   «  fa\2        2 

where  a  »  (•<- 1)~  •     As  l°ng  as  H     is  different  from  zero,  m  is  a  simple  zero  of 
det(C).     The  solution  corresponding  to  such  an  m  is  a  contact  discontinuity.     The 
nature  of  this  solution  will  be  discussed  later  in  subsection  2.U.     The   solutions 
corresponding  to  the  zeros  of  the  factor  (m  X    -  nH^)  are  called  transverse  shocks   , 
de  Hoffman  and  Teller  call  such  shocks  symmetrical  shocks. 
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and  the  solutions  corresponding  to  the  larger  and  smaller  zeros  of  the  remaining 

f—  \2         2 
factor,   assuming  (  H  )    -  H     does  not  vanish,   are  called  fast  and  slow  shocks  re- 
spectively.    Among     the  various  properties  of  these   shocks  the  following  are  use- 
ful for  our  purposes.     First,   fast  and  slow  shocks  are  compressive,  that  is,   the 
density  and  pressure  behind  the  shock  exceed  the  corresponding  quantities  ahead 
of  the   shock,  while  in  a  transverse   shock  only  the  transverse  magnetic  field  and 
transverse  velocity  suffer  a  discontinuity.     Second,   the  magnetic  field  strength 
|H.     j      rises  across  a  fast  shock,  drops  across  a  slow  shock,  and  remains  constant 
across  a  transverse  shock   .     Third,  under  certain  circumstances   switch-on  shocks 
exist.     Switch-on  shocks  are   shocks  having  a  vanishing  transverse  magnetic  field 
component  in  front  but  nevertheless  a  non-vanishing  transverse  magnetic  field 
component  behind  the  shock.     Clearly  switch-on  shocks  must  be  fast  shocks.     In 
Section  3,  where  we   shall  treat  switch-on  shocks  in  greater  detail,  we  shall 

prove  that  the  normal  flow  velocity,  u  -U  ,  is  equal  to  the  so-called  Alfven 

2   -1  1/2 
speed  b  «(uH  p     )         behind  the  shock.   It  should  be  mentioned,   in  this  connection, 

that  when  the  transverse  magnetic  field  vanishes  both  in  front  and  behind,  the  dis- 
continuity surface  the  relations  B,    -  B,    reduce  to  the  well-known  Rankine-Hugoniot 
relations  for  gas  dynamical  flow  with  a  constant  transverse   flow.     The  velocity 
of  the  gas  shock  in  this  system  is  in  no  way  related  to  the   Alfven  speed  behind 
the   shock.     It  is  nevertheless  useful  to  divide   gas  shocks  into  two  classes 

according  as  the  normal  flow    velocity   (u     ,-U    )  behind  these  shocks  is  equal  to 
e  iifl     n 

or  greater  than  the  Alfven  speed  behind  the  shock,  or  less  than  this   speed.     In 
the  former  case  the  gas  shocks  are  called  fast  gas  shocks;   in  the  latter  case 
they  are  called  the  slow  gas  shocks. 

2.3       Hydromagnetic  simple  waves 

In  a  given  flow  region  R  of  the   (x,t)-plane,  simple  wave   solutions  of 

equation  (2.1)  are  continuous,  continuously  different iable   solutions  of  a  special 

_       .      ___ 

Since  H.   is  discontinuous  across  a  transverse  shock  it  follows  that  the 
tr 

direction  of  ff   must  change. 
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type.     They  may  be  characterized  by  the  requirement  that  H,  u,  P,  p,  and  S  be 

constant  on  a  one-parameter  family  of  lines  which  cover  R.     This  family  of  lines 

may  be  written  in  the  form 

x     -    Ut  +  £,  U     -     B(s), 

if  not  all  the  lines  pass  through  the  origin  of  the  (x,t)-plane,  or  in  the  form 

x  -  Ut     >     U  -  U(£) 

if  they  do.  A  simple  wave  whose  flow  region  is  covered  by  the  latter  family  are 
referred  to  as  centered  simple  waves.  In  a  centered  simple  wave  S",  u,  p,  p  and 
S  are  necessarily  functions  of  U(»  x/t)  alonei  otherwise  these  quantities  are 
necessarily  functions  of  the  parameter  £(■  x-Ut).  It  can  be  shown  that  a  flow 
region  adjacent  to  a  region  of  constant  state  is  always  a  simple  wave  and  that 
the  transition  from  a  simple  wave  region  to  a  region  of  constant  state  takes  place 
across  a  line  of  the  type  discussed  above.  It  is  for  these  reasons  that  simple 
waves  are  useful  building  blocks  in  constructing  solutions  to  flow  problems. 
For  solutions  of  simple  wave  type,  (2.1)  reduces,  on  setting 

(2.12)  U  -  ux  *  c, 

to  the  following  system  of  ordinary  differential  equations: 
\j  *cdHy*Hydux-Hjcduy     -     0 

D.,  TcdH+Hdu-Hdu       -0 

1,Z  Z  Z      X  X      z 

E^  +  pcdux  +  dp  ♦  j.  nd(f2)     -    0 

(2.13) 

Dg  *  pcduy  -  ui^  dHy     -     0 

D2,z  *  pCdUz  "  ^z     "     ° 

D-  ♦  cdp     +     pdu       ■     0 

3  ^ 

D.  *  cdS     -    0, 

*A  state  of  flow  is  referred  to  as  constant  in  a  region  of  the   (x,t)-plane  if  all 
quantities,  p,  p,  H,  u,  S  are  constant  throughout  the  region. 
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where  the  differentials  involving  the  parameter     £  have  been  multiplied  out. 
For  non-trivial  solutions  the   quantity  c  must  be  chosen  in  such  a  manner  that 
the  determinant  of  the  system 

(2.1U)  det(D)   -  ±  pc2(pc2-ui£)r(pc2  -  nH2)(c2-a2)-  c2u-(H2  ♦  H2)l 

-1  1/2 
vanishes.     Here,  a  »   (yPP     )         is  the  sound  speed.     The   zeros  of  this  deter- 

minantal     expression  are  referred  to  as  the  characteristic  disturbance  speeds. 
The  zero  c  -   (\±kTp~   )         *  b  is  called  the  Alfven,  or  the  transverse   disturbance 
speed,   the   smaller  zero,  c   ,  of  the  factor  in  brackets  is  called  the  slow  dis- 
turbance   speed,   and  the  larger  zero,   cf,  the  fast  disturbance  speed.     The  velo- 
cities U  -  u     ±  b  are  accordingly  called  Alfve'n  (or  transverse)  characteristic 
velocities,  and  U  ■  u     ±  c   ,  U  -  u     ±c_  are  called  the  slow  and  fast  character- 

X  S  A  1 

istic  velocities.  According  to  which  speed  c  enters  into  the  solution  of  (2.13), 
the  simple  wave  is  transverse,  slow  or  fast.  The  wave  corresponding  to  c  ■  0 
is  called  a  contact  layer. 

A  solution  of  the  system  (2.13)  corresponding  to  the  upper  (or  lower)  sign 
and  hence  to  U  ♦  c  (or  U  -  c)  is  referred  to  as  a  forward-facing  (backward- 

X  A 

facing)  simple  wave.     As  is  the  case  with  forward-facing  (backward-facing)  shocks, 
fluid  particles  enter  the   simple  wave  region  from  the  right  (left)  if  the  simple 
wave  is  forward-facing  (backward -facing).     This  remark  is  a  consequence  of  the 
fact  that  u     +  c   (or  u     -  c)  is  larger   (smaller)   than  u ,  the  longitudinal  velo- 

XX  X 

city  of  the  fluid. 

In  the  following  survey  of  the  properties  of  the  various  simple  wave  solutions 
it  will  be  assumed  that  H  f   0.  Across  a  contact  layer  only  p  and  S  may  vary, 


the  other  quantities  being  necessarily  constant.  Across  a  transverse  wave,  p,  p, 
S  and  u  are  constant.  The  transverse  component  of  the  magnetic  field  remains 
constant  in  magnitude  but  may  change  its  direction  arbitrarily  but  continuously 
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with  xj  the  transverse  component  of  the  velocity  u.  is  given  simply  by 

u.  ■  a  ±  cH"  M.  where  a  is  a  constant  transverse  vector.  Across  a  fast 
tr    *   x  tr> 

or  slow  simple  wave  S  is  constant j  the  magnitude  of  all  other  quantities  may 

vary.  A  detailed  description  of  these  waves  will  be  given  in  Section  U.  It  is, 

however,  worthwhile  noting  here  that  the  direction  of  H.  remains  constant  across 

fast  and  slow  simple  waves  and  that  du.  is  parallel  to  H.  ,  i.e.,  du,  x  H. 

vanishes.  The  fact  that  du".  x  H   vanishes  is  a  consequence  of  IL   and  D,  . 

"tr    tr  c , y     c  jZ 

The  validity  of  the  first  statement  is  indicated  by  the  following:  Multicly 

equation  D,   by  H  and  equation  T)0       by  H  and  subtract  the  former  equation 
!.,y    z  £,z    y 

from  the  latter.  Then  eliminate  du  and  du  from  the  result  by  means  of  D0 

y      z  2,y 

and  Do  •  We  then  havs 
2,z 

(c  -  uH^pc)"1)  (HdH  -  H  dH  )  -  0. 
x  z  y    y   a 

Since  it  can  be  shown   (see  Section  il]  that  the  first  factor  vanishes  at  only  one 

point  of  a  slow  or  fast  simple  wave,  we  conclude  that  the  H  ■  const *H.  It  is 

assumed  here  that  neither  H  nor  H  vanish  identically  in  the  fast  or  slow  simple 

y  z 

wave  region  under  consideration.     Corresponding  relations  exist  for  shock  waves. 
In  fact  it  is  easily  verified  that  the  direction  of  H.     behind  a  shock  is  collinear 
with  that  ahead  provided    |ff     |   does  not  vanish  ahead  of  the   shock    (cf.   (3.1)]. 
Moreover,   since  the  transverse  part  of  Cp  is 

■Bi,]   -MHn   Btrl     "     ° 

it  follows  that    fu      "1  x   fET     j  vanishes.     In  particular,  independently  of  whether 

H.     vanishes  ahead  of  the   shock  or  not,  if  the  transverse  magnetic  field  is  assumed 
tr 

to  be  parallel  to  the  y-axis  we  may  conclude  that  u     is  constant  across  a  shock. 
Similarly,  if  B.     is  parallel  to  the  y-axis  then  u    is  constant  across  fast  and 
slow  simple  waves. 
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2.U     Formulation  of  the  problem 

We  turn  now  to  the   formulation  of  the  problem.     In  the  following  it 
will  be  assumed  that  H    >  0,  unless  otherwise  indicated.     Since  we  are   dealing 
with  one-dimensional  motion,  discontinuity  surfaces  are  necessarily  planes 
normal  to  the  x-axis.     Let  L  be  an  initially  stationary  discontinuity  surface 
located  at  x  »  0  across  which  there  exists  at  time  t  ■  0  a  shear  flow  discon- 
tinuity.    By  a  shear  flow  discontinuity  we  mean  one  in  which  the  transverse 

component  of  the  velocity  u       jumps  without -changing  direction,  the  variables 

or 

u  >   E>  P»  P>  S  remaining  continuous  and  constant  for  all  x.  In  particular  i« 
is  assumed  that  at  time  t  ■  0 


(2.15) 


uy  "  uy,o  >  °'     x  >  °>        uy  "  "uy,o>      x  <  °  ' 


P  -  P0«   P  -  P0,   ux  -  uz  -  0,   Htr  -  0,   S  -  S0 


[see  Figure  l].     It  should  be  observed  that  there  is  no  loss  in  generality  in 
assuming  that  u     and  u     vanish  and  that  the  transverse  speeds  on  either  side  of 
the  discontinuity  surface  are  equal,  since  a  more  general  prescription  of  the 
initial  motion  may  be  reduced  to  this  case  by  means  of  a  galilean  transform ationj 
it  is  easily  verified  that  the  equations  of  motion  (A)  and  the  shock  conditions 
(B)  are  invariant  under  such  transformations. 

Since  the  normal  flow  velocity  u     vanishes  initially  and  since  L  is  assumed 
to  be  stationary  initially,  it  follows  that  the  flux  through  L,  m,  vanishes.     The 
shear  flow  discontinuity  therefore  appears  to  be  some  form  of  hydromagnetic  contact 
discontinuity.     It  will  now  be  shown  that  such  is  not  the  case  -  more  explicitly, 
that  a  shear  flow  discontinuity  is  not  an  admissible  hydromagnetic  motion  if  the 
longitudinal  magnetic  field  H     is  different  from  zero.     To  this  end  let  ue  set  m 
equal  to  zero  in  the  Rsnfcine-Hugoniot  conditions,   (2.U).     Equations  (2.h)-B     and 
B.  and  the  entropy  condition    then  yield  the  following  relations t 
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(2.16)  pj  -  0,  [un]  -  0,        [p]  >  0,  [S]  >0. 

These  relations  in  conjunction  with  (2.U)-Bp  and  B.    then  imply 

(2.i7)  [p]    -   o,        prtP]  -  0. 

Finally,  the  fact  that    [h.    J  and    [h  ]  vanish  implies,   in  virtue  of  (2.U)-B,,   that 

(2a8)  Hn   [fftJ    =     °    ' 

where  u       is  the  transverse  component  of  the  fluid  velocity  u. 

Now  the   shear  flow  discontinuity  (2.15)   clearly  satisfies  the  corditions 
listed  in  (2.16)  and   (2.17).     Equation  (2.18)   is,   however,  net  satisfied  since 
we  have  assumed  that  H     *  H    >  0.     The  shear  flow  discontinuity  in  question  is 
therefore  not  an  admissible  hydromagnetic  motion. 

Note  that  when  H     is  permitted  to  vanish  u.      may  be  discontinuous  across 
L;  in  this  case  the  problem  reduces  to  the  one -dimensional  gas  dynamical  prob- 
lem of  a  shear  flow  discontinuity  in  a  homogeneous  medium.     A  shear  flow  dis- 
continuity is,  under  these  circumstances,   an  admissible  gas  dynamical  motion 
and  may  persist  for  all  time.     On  the  other  hand,  when  H    j  0.  the  initial  dis- 
continuity must  be  resolved  by  a  suitable  motion  of  the  medium. 

To  gain  some  insight  into  the  cause  of  the  motion  consider  instead  of  a 
discontinuity  a  thin  shear  flow  layer  across  which  the  transverse  component  of 
velocity,  u     say,  varies  smoothly  from  its  value  u     ■  -u     Q  o/>  the  left-hand 
side  of  the  layer  to  the  value  u     «  u     Q  on  the  right-hand  side.     To  simplify 
matters  assume  that  u     vanishes  everywhere.     In  addition  assume  that  H     ■  H„  ■  0 
initially,  and  H    >  0.     For  the  sake  of  concreteness  let  us  suppose  that  u     is 

*  x  y 

an  odd  function  of  x  throughout  the  layer  which  we  shall  assume  to  be  symmetrically 
disposed  with  respect  to  the  (y,z)-plane.  From  (2.2)  we  conclude  that  the  electric 
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field   ,E-E     z     ■  -(ua  H  )*\  .is  an  odd  function  of  x  which  varies  smoothly 
'  %     o  Vl    y  x     o »  * 

from  (till     ~H  )z     on  the  left  face  of  the   slab  to  {-vox     nE   )z     on  the  right 
y,0  x     o  y>0  x'  o 

face  of  the  slab.     The  existence  of  this  initial  electric  field  results  in  the 
growth  of  a  transverse  component  of  the  magnetic  field  in  the  layer.     For  t 
sufficiently  small  this  transverse  magnetic  field,   FL    ,   is  given  by    see   (2.2 )-A 

3u 

H.        «     H     J    ty   . 
tr  x  ax  *  o 

Thus  H       is  directed  only  along  the  positive  y-axis  and  is  an  even  function  of 
x,  since  u    is  odd.     The  current 

3H  2 

7    -     *  x  fftr     "     1ST    *o  "  <H*  7T  *>»o 

oX 

is  then  an  odd  function  of  x  which  is  negative  when  x  >  0  and  is  directed  along 
the  z-axis.  The  force  per  unit  volume, F,  is  then  given  by  the  expression 

X  3X   3x2      ° 

This  force  is  directed  along  the  positive  x-axis  when  x  >  0  and  along  the 

negative  x-axis  when  x  <  Oj  thus  it  tends  to  produce   separation.     Note  that 

we  could  have  assumed  that  the  longitudinal  magnetic  field  was  directed  along 

the  negative  x-axis  or  that  the  initial  transverse  velocity  was  u   (-x)  without 

changing  the  character  of  F. 

One  can  readily  infer  from  the  above  discussion  that  the  greater  u    n,   (H 

remaining  fixed)  the  greater  the  force  tending  to  produce   separation.     If  u     . 

is  made  large  enough  one  might  expect  the  layers  to  separate  altogether,  leaving 

a  vacuum  in  between.      That  this  situation  may  actually  arise  will  be  seen  later. 
In  determining  the  motion  which  resolves  the  initial  discontinuity  we  admit 

solutions  of  the  following  type.     First,  we  admit  only  those   solutions  which  may 

be  constructed  from  combinations  of   shocks,   simple  waves  and  constant  states. 

#  _   ~    — — ■   —  — ^-——^-—— —————— ——_—_— _______ — —  ■ 

The  vectors  xq,  y  and  z  constitute  the  unit  orthogonal  tried, 
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Second,  admissible  solutions  are  required  to  depend  continuously  on  the  para- 
meters u  _,  H  ,  p  and  p  -  that  is.  on  the  parameters  which  characterize  the 
y,0   x   o     o        ' 

initial  state.  Third,  we  require  that  the  fluid  velocity,  u,  be  an  odd  function 
of  x  while  H,  p  and  p  are  even  functions  of  x.  Fourth,  we  shall  from  the  outset 
limit  the  class  of  solutions  to  those  for  which  the  transverse  component  of  velo- 
city is  parallel  to  the  y-axis,  i.e.,  to  those  for  which 

(2.19)  u   -  0,     -co  <  x  <  oo,     t  >  0. 

As  already  mentioned  in  the  introductory  section,  the  resolution  of  the 
initial  discontinuity  may  be  effected  by  two  fast  shocks  (switch-on  or  fast  gas 
shocks  depending  on  the  thermodynamic  state  of  the  undisturbed  medium  and  the 
magnitude  of  the  initial  discontinuity)  followed  by  slow  centered  simple  waves. 
That  the  shocks  in  the  resolution  must  be  switch-on  or  pure  gas  shocks  follows 
from  the  fact  that  the  transverse  magnetic  field  is  zero  initially  and  hence  is 
absent  in  front  of  the  advancing  shocks.  The  choice  of  slow  simple  waves  is  re- 
lated to  the  requirement  that  the  shocks  be  switch-on  or  fast  gas  dynamical  shocks. 
It  will  be  recalled  that  the  speed  |U  |  of  both  types  of  shocks  are  at  least 
equal  to  the  Alfven  characteristic  speeds  |u  *  b  |  behind  the  shocks  (see  the 
end  of  subsection  2,2]»  In  slow  (or  fast)  simple  waves  the  characteristic 
speeds  |u  ±  c|  are  at  most  (or  at  least)  equal  to  the  corresponding  Alfven 
characteristic  speeds, as  we  shall  show  in  Section  U.  Of  the  two  types  of  simple 
waves  clearly  only  the  slow  simple  waves  "  follow"  switch-on  or  fast  gas  shocks. 

The  third  and  fourth  requirements  on  our  class  of  solutions  are  easily  seen 
to  be  compatible  with  the  initial  conditions,  the  continuum  equations  and  the 
shock  conditions.  Moreover,  they  are  consistent  with  the  heuristic  results 
obtained  in  our  analysis  of  the  shear  flow  layer.  Note  that  the  vanishing  of 
u  everywhere  implies  that 

(2.20)  H  -  0,    -oo  <  x  <  oo,   t  >  0 
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in  the  fast  shock  slow  simple  wave  resolution.  This  assertion  follows  from  the 

relation,  du"   x  ff   -  0,  which  holds  in  slow  and  fast  simple  waves  and  from  the 

'       tr         tr 

relation,    [u     ]  x    [X..J   "  °j  which  holds  across  any  admissible  discontinuity. 

The  converse  is  also  true  since  u  ,  by  hypothesis,   vanishes  in  front  of  the  ad- 

z 

vancing  shocks.  Thus,  under  the  present  circumstances,  the  vanishing  of  Hz 

everywhere  is  equivalent  to  the  vanishing  of  u  everywhere.  It  should  be  recalled 

z 

here  that  the  direction  of  ff      in  slow  (and  fast)   simple  waves  is  constant  across 

tr 

the  wave.     The  requirement  that  H.      be  parallel  to  the  y-axis  in  such  a  wave  is 
therefore  not  a  very  severe  restriction. 

Consider  now  the  region  between  the  slow  simple  waves  in  the  resolution  of 
the   shear  flow  discontinuity  discussed  above.     This  region,  being  bordered  by 
simple  wave  regions,   is  clearly  one  of  constant  state.     Since  u  is  odd  we  must 
have 

(2.2i)  u-0,  u-0,  t>0 

v  x         '  y         ' 

unless  p  vanishes  throughout  the  region.  One  must  admit  the  alternative  p  ■  0 
in  order  to  take  into  account  the  possibility  of  cavitation.  As  indicated  above, 
the  likelihood  that  cavitation  occurs  increases  with  the  magnitude  of  the  initial 
discontinuity,  assuming  the  remaining  parameters  in  the  undisturbed  medium  are 
fixed.  When  cavitation  occurs  it  is  necessary  to  replace  conditions  (2.21)  by 
the  conditions 

p     -     0 


(2.22) 

Ly~x      ~x 


u_H„  -  u_H     *  0 


at  the  endpoints  of  the  vacuum  zone. 


The  second  equation  simply  states  that  the  electric  field, 

(2.23)  E  -  Mff  x  u     -    n£yHx  -  «J  zQ  , 

vanishes  at  the  endpoints  of  the  vacuum  zone.  This  result  is  a  consequence  of 
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the  fact  that  E  vanishes  throughout  the  vacuum  zone,  as  we  shall  now  show* 
Maxwell* s  equations  in  the  vacuum  zone  are 

3E.  (x,t)      m     . x  3E  (x,t) 


xo    * 

SEtJ((x,t) 

ax 

*0  * 

3Htr(x,t) 

dx 

3H 
(2.25)     x„  x   -S^ 0  ,  ar  "  °» 


where  as  usual  we  have  assumed  that  the  displacement  current  may  be  neglected  and 

that  all  variable  quantities  depend  on  x  and  t  alone.  In  accordance  with  our 

previous  notation,  E~.  and  ff.  are  the  transverse  electric  and  magnetic  field 

x»r  x»r 

components,  E     and  H.  the  longitudinal  components  of  the  electric  and  magnetic 

field,  ff  the  total  magnetic  field,  and  x    the  unit  vector  directed  along  the 

positive  x-axis. 

The  two  equations  in  (2.25)  imply  that  H  is  independent  of  x  and  t  and  that 

ff.  is  independent  of  x.  From  (2.2U)  we  then  conclude  that 
tr 

(2.26)  Ex  "  kl' 

3Htr 

(2.27)  E"tr-  nxox  -5T-  X+*2, 

where  k1  is  a  scalar  quantity  and  k?  a  transverse  vector  quantity  which  depend 
on  t  alone.     Now  it  can  be  shown    [see  Appendix  I]   that  the  Rankine-Hugoniot 
conditions  at  the  boundary  of  a  perfectly  conducting,  electrically  neutral, 
compressible  fluid  which  expands  into  a  vacuum  under  action  of  magnetic  and 
pressure  forces  are*  assuming  H  •  n  f  0, 

(2.28)  [p]   -  0,  [If]     -     0,  [ff]      -     0  . 

The  above  relations  imply  that  the  density  is  continuous  at  this  boundary  and 
vanishes  there  and  that  surface  current  layers  and  charge  layers  (single  or 
double)  are  absent.  We  observe,  at  this  point,  the  fact  that  the  boundary  of 
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a  simple  wave  region  consists  of  linear  characteristics;  in  the  present  case, 

since  the  slow  simple  waves  are  centered,  these  linear  characteristics  pass 

through  the  origin  of  the  (x,t)-plane.  In  addition,  we  recall  that  the  state 

variables,  and  ff*.   in  particular,  are  constant  along  such  characteristics.  Thus 
tr 

the  direction  and  magnitude  of  H  are  constant  at  the  moving  boundary  of  the 
simple  waves.  Since  H   is  continuous  at  the  boundaries  separating  the  vacuum 
zone  from  the  simple  wave  zone  we  conclude  that  H,  ,  in  addition  to  being  in- 
dependent of  x  in  the  vacuum  zone,  is  also  independent  of  t.  It  follows  immedi- 
ately from  (2.27)  that 

(2.29)  Ftr  •  K2  . 

But  u  is  odd  and  ff  is  even.     Therefore  if  E  -  nH  *u  is  the  electric  field 

at  one  endpoint  of  the  vacuum  zone,  E*    -  -uHvuis  the  electric  field  at 

the  other.     Employing  this  fact,   the  fact  that    [  E  ]  vanishes  at  the  boundaries 

of  the  vacuum  zone,  and   (2.29),  we  conclude  that 

(2.30)  t"tr     -     0. 

In  addition,  we  note  that  the  electric  field  in  the  simple  wave  region  is  entirely 
transverse  since  both  u  and  H  lie  in  the  (x,y)-plane.  It  follows  from  (2.26) 
that  E  vanishes.  Thus  we  have  shown  that  the  electric  field  vanishes  throughout 
the  vacuum  zone. 

It  should  be  observed  that  the  condition  u  ■  0  imposed  by  Friedrichs  [see 
[7]  pp.  U8-U9]  on  the  transverse  component  of  the  fluid  velocity  at  the  linear 
characteristics  which  bound  the  vacuum  zone  -  that  is,  at  the  so-called  "tails" 
of  the  simple  waves  -  is  incorrect  and  must  be  replaced  by  the  less  restrictive 
condition  given  in  (2.27).  This  change  affects  the  qualitative  picture  of  the 
motion  given  by  Friedrichs  only  when  the  initial  discontinuity  2u  n     is  large 
enough  to  produce  cavitation, and  then  only  in  the  following  manner.  The  fluid 
velocity  at  the  tails  of  the  simple  waves  has  transverse  components  in  addition 
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to  the  escape  velocities  in  the  longitudinal  direction.  At  each  boundary 
the  second  equation  of  (2.22)  governs  the  relation  of  the  transverse  velocity 
components  to  the  longitudinal   escape  velocities. 

It  should  be  pointed  out  that  this  change  also  affects  the  relation  between 
the  transverse  magnetic  field  developed  behind  the  waves  and  the  magnitude  of  the 
initial  discontinuity,  2u,  _,  when  u  n   is  large.  Equations  (13.5)  and  (13.7) 
of  [7]  must  be  replaced  by  (1.2)  of  the  present  report  (see  also  (5.37)). 

Some  remarks  are  in  order  here  regarding  the  choice  of  pairs  of  waves  in 
the  resolution  of  the  initial  discontinuity.  The  choice  of  pairs  of  waves  is 
motivated  by  the  following  considerations.  One  expects  shock  waves  to  be  com- 
pletely determined  by  the  prescription  of  one  quantity,  the  pressure  say, 
behind  the  wave,  provided  the  state  in  front  of  the  wave  is  known.  The  same 
may  be  said  of  simple  waves.  One  may  therefore  hope  to  be  able  to  prescribe 
a  pair  of  such  quantities  behind  a  combination  of  a  simple  wave  and  a  shock. 
In  particular,  one  may  hope  to  adjust  these  parameters  in  such  a  way  that  (2.21) 
or  (2.22)  is  satisfied.  This  procedure  is  carried  out  in  detail  in  Section  5 
where  the  dependence  of  these  parameters  on  the  initial  discontinuity,  2u  Q, 
and  on  the  thermodynamic  state  is  determined. 

It  should  be  understood  from  the  outset  that  no  claims  are  made  with  regard 
to  the  uniqueness  of  our  solution.  The  solution  does,  however,  appear  to  be 
the  only  ons  satisfying  the  initial  conditions  if  we  limit  the  class  of  admissible 
solutions  to  that  discussed  above  and,  moreover,  impose  the  additional  requirement 
that  we  employ  at  most  two  pairs  of  shocks  or  simple  waves  or  a  pair  of  shocks 
and  a  pair  of  simple  wave a.  The  validity  of  this  statement  may  be  tested  by 
attempting  to  construct  from  the  basic  solutions  discussed  in  subsections  2.2 
and  2.3  an  alternative  solution  satisfying  all  the  necessary  requirements. 
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3.  Switch-on  and  fast  gas  dynamical  shocks 
3.1  Introductory  remarks 

As  already  mentioned,  svitch-on  and  pure  gas  dynamical  shocks  play  a  role 
in  the  resolution  of  the  initial  shear  flow  discontinuity.     Switch-on  shocks,  it 
will  be  recalled,   are  fast  shocks  which  have  a  non-vanishing  transverse  component 
of  the  magnetic  field  behind  the  shock  in  spite  of  the  absence  of  such  a  compon- 
ent in  front  of  the  shock.     Fast  gas  shocks  are  pure  gas  dynamical  shocks  satis- 
fying the  condition  that  the  normal  flow  velocity  behind  the   shock  (u     ,   -  U")  is 

u«  JL  n 

greater  or  equal  to  the  Alfve'n  disturbance   speed,  b,    ■  (uJT  p7   )       ,  behind  the 

shock.     This  section  is  devoted  to  discussing  the  properties  of  these  shocks. 

We  shall  show  that  svitch-on  shocks  exist  provided    s     (  -  vp  /H  ) .  the  sauare 

o     iro  x  '     n 

ratio  of  the  sound  speed  a   -  (yp  p"  )  '    \   and  the  Alfve'n  disturbance  speed 
b   •  ((jfl  p"  )  '    in  front  of  the  shocks,  is  less  than  unity  and  provided  the 
shock  strengths  measured  in  terms  of  the  pressure  ratio,  is  below  a  certain  criti- 
cal value.  When  these  conditions  are  satisfied  it  will  be  shown  that  the  state 
behind  the  switch-on  shock  is  known  once  the  state  in  front,  the  direction  of  the 
transverse  magnetic  field,  and  the  pressure  ratio  is  prescribed.  That  the  state 
behind  a  gas  shock  is  determined  by  the  state  in  front  and  the  pressure  ratio  is 
a  well-known  factj  we  shall  therefore  merely  quote  the  results  in  a  form  suitable 
for  subsequent  purposes.  When  s  <  1  we  shall  prove  that  the  gas  shock  is  fast 
if  and  only  if  the  pressure  ratio  is  above  the  critical  value  alluded  to  above. 
In  addition,  we  shall  show  that  as  the  pressure  ratio  passes  through  the  critical 
value,  the  state  behind  a  switch-on  shock  having  a  given  state  in  front  passes 
continuously  into  that  of  the  fast  gas  shock  having  the  same  state  in  front. 
It  might  be  worthwhile  adding  here  that  only  after  these  shocks  are 
combined  with  the  appropriate  simple  waves  in  forming  the  solution  of  the  shear 
flow  problem  is  it  possible  to  fix  the  direction  of  the  transverse  magnetic  field 
behind  the  switch-on  shocks  and  to  relate  the  shock  strength  to  the  strength  2u 
of  the  initial  diacontinuity • 
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3.2     Switch-on  shocks 

The  analysis  of  switch-on  shocks  is  based  on  the  following  jump 
relations: 


(3.1) 


f(»2r    -  hhJ;)  irtrJ  -    0,  m   f    0 


(3.2)  m2+   £  +  ^(Htr)2]    IXT1     "     ° 

(3.3)  (*2t    -  MB*)  J  £  t     -  Oj.  p  -  njSf2  ♦  n(f)J 


&     J(«2-H2U     0. 


The  first  relation  is  obtained  from  the  transverse  part  of  B,,  namely 

(3.11)  £*    fftr-Hnv]  "    °* 

on  substituting  for    u      I  from  the  transverse  part  of  the  momentum-conservation 
shock  condition  B«,  that  is, 

(3.5)  f"  V  -  ^  g]     '    0- 

Equation  (3.2)  follows  from  the  longitudinal  part  of  B2  and  the  flux-conservation 

jump  condition,  B.,.     The  last  relation,  is  the  equation  which  results  from  equating 

the  last  factor  in  the  right-hand  side  of  (2.11)  to  zero.     In  accordance  with  the 

conventions  of  the  previous  section  0     and  Q.   are  the  values  of  any  quantity  Q 

ahead  of  and  behind  the  switch-on  shock  3    -     (Q,   +  Qg)/2  and   [Q]   »  Q,   -  QQ. 

Since  H".     _  vanishes  in  front  of  a  switch-on  shock  we  have  from  (3.1) 
tr,u 


(3.6)  (m2t  x  -  hhJ)  WtrA    -    0. 
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A  necessary  condition  for  the  existence  of  a  switch-on  shock  is  therefore 

(3.7)  «2f  i  -  MHj  -  0. 

2 

Eliminating  m  from  (3.2)  by  means  of  this  relation  we  have 

(3.8)  2T(fftr,l)2  "  -M  *  ^pl  M. 
Since  for  switch-on  shocks 

(3.9) 

M"  [fftr]2  ■  «*«!>*. 

then  on  eliminating  (3",   ,  )  from  (3.3)  by  means  of  (3.8)  ve  obtain,  employing 
(3.7)>  an  equation  involving  only  pQ,  pQ  and  the  ratios  (p±-P0)/Pq»  (P1-PQ)/P0 
This  equation  when  solved  for  the  ratio  (p-,-p0)/p  yields 

or  alternatively,  when  solved  for  the  ratio  (p,-p  )/p  , 


It  is  apparent  from  the  fact  that   (p,-p  )/p     is  quadratic  in  (p1-p  )/p    that  the 
sign  before  the  square  root  terra  in  (3.10)  could  also  be  negative.     However,   such 
a  choice  of  sign  would  lead  to  negative  values  for  the  ratio   (p,-p  )/p  ,  that 
is,  to  switch-on  shocks  which  are  not  compressive.     Since  all  shocks,  with  the 
exception  of  transverse  shocks  are  compressive,  the  choice  of  the  positive  sign 
is  required. 
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Substituting  the  expression  for   [p]   in  (3.IJ)  into  (3.8)  we  find 


(3.12) 


<±<AS  •>*&{- &&&*%]} 


The  fact  that  the  right-hand  side  of  this  equation  is  positive  for  a  switch-on 
shock  results  in  the  following  restriction  on  the  range  of  variation  of  the  ratio 

(pl-po)/po: 


(3.13)  0    <     ?P><    £     [l-(YP0M)] 


This  relation,  in  turn,  implies,  in  virtue  of  (3.10) 

uH2 

(3.iu)  o   <   -^  <   Sj.   (  -5  )  |i  -  (yp0/mh2)]  . 

o 

Within  this  range,  the   quantity  (H.    /H  )     increases  from  zero  to  the  maximum 

(3.15)  (VAlL    ■    ^T    ^-(rPo/^)]2 

which  occurs  when 

P-.-P. 


(3.16)         (Xo).  Jj-ji.C^/^)] 


and  then  decreases  to  zero  as  the  critical  value  of  the  excess  pressure  ratio 


w 


(3.XT)  (^)crit    ■    Sc 


&S  E-^M'] 


o  '  ,jro 

is  approached     cf.  (3.1U)     . 

We  have  thus  far  shown  that  the  relations  (3.6)  and  (3.13)  or  (3.U»)  are 
necessary  conditions  for  the  existence  of  switch-on  shocks.     In  addition  to  these 
conditions  we  note  that  the  braced  term  in  (3.13)  or  (3. Hi)  is  positive  and  that 
hence 
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YP 
(3.18)  °    <     i. 

x 


Conversely,  if  a  shcck  has  a  vanishing  transverse  magnetic  field  ahead  and  is 
such  that  the  relations  (3. 7), (3.13)  or  (3.1U)  and  (3.18)  are  satisfied,  then 
a  non-vanishing  transverse  magnetic  field  exists  behind  the  shocks. 

Consider  now  the  problem  of  determining  the   relationship  between  the  state 
(1)  behind  a  switch-on  shock  and  the  state    (0)  in  front.     Thus  far,   states  have 
been  described  by  the  variables  p,  p,  u~  and  H,    ,  it  being  assumed  that  H     has 
been  specified  once  and  for  all.     By  way  of  preparation  for  subsequent  consider- 
ations, however,   it  is  convenient  to  introduce  the  alternative  set  p,  s,  b,  u 

and  ff.    ,  where, 

tr*  ' 

(3.19)  b     -     IP  b     -     (nH*/p)1/2     . 

x 

Described  in  physical  terms  b  is  the  Alfven  speed  and  s  is  the  square  ratio  of 


the  sound  speed    y  yP0~     and  the  Alfven  speed. 

In  the  following  it  will  be  assumed  that  we  know  which  region  in  the 
(x,t)-plane  is  the  region  behind  the  shock     and  employ  the  ratio   (s.-s   )/s 
("  (p-i-P  /p   ))  as  a  measure  of  the   shock  strength.     It  will  be  recalled    (cf. 
Section  2.2}  that  the  normal  to  the  shock  front,  n,  is  directed  into  the  region 
behind  the  shock j  n  is  therefore  directed  along  the  negative  or  positive  x-axis 
according  as  the   shock  is  forward  or  backward  facing. 

It  will  now  be   shown  that  both  the  state   (1)  behind  the  shock  and  the 
velocity  U     of  the   shock  front  are  completely  determined  whenever  the   state   (0) 
ahead  of  the  shock,  the  shock  strength  (s,-s  )/s  ,  and  the  direction  of  the  trans- 
verse magnetic  field  are  specified.     The  necessity  for  prescribing  the  direction 

Since  switch-on  shocks  are  compressive  it  is  only  necessary  to  specify  which 
region  of  the   (x,t)-plane  is  the  region  of  higher  pressure. 
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of  the  magnetic  field  is  indicated  by  the   following  argument.     Consider  first 
the  general  situation  where  the  transverse  magnetic  field  does  not  vanish  in 
front  of  the  shock.     The  hydromagnetic  Rankine-Hugoniot  conditions  B,,...,B. 
represent  seven  conditions  on  the   fifteen  variables  p,,  p,,  H.     ,,  u.      ,,  u     ,, 
p  ,  p  ,  ff.      0,  u.     0,  u    Q,  U   .     When  the  state  in  front  is  prescribed,  the 
total  number  of  conditions  is  fourteen;  these  conditions  yield,  in  general,  a 
one-parameter  family  of  states  behind  the  shock.     Assume  now  that  the  shock  is 
a  switch-on  shock.     In  this  case  the  four  Rankine-Hugoniot  conditions,   B,   and 
the  transverse  part  of  B2   [cf.   (3.U)  and  (3.5)3*  degenerate  into  three  condi- 
tions,  namely  those  given  by  the  relations   (3.5)  and   (3*7).     If  we  prescribe 
the  state  in  front  we  obtain  thirteen  conditions  on  the  fifteen  variables 
mentioned  above.     To  determine  the  state  behind  the  shock  it  is  necessary  to 
specify,  in  addition  to  the   shock  strength,  another  parameter.     We  may  choose 
this  parameter  to  be  the  direction  of  the  transverse  magnetic  field  behind  the 
shock. 

Without  loss  of  generality  we  may  assume  that  H.      is  parallel  to  the 
y-axi8.     The   state   (1)  and  the  velocity  of  the  shock  front  U    are  then  given 


by  the  following  relations: 

i  p0  y=T   )         L  y  so       so 


E2 


^/-^-[<^>-oj 


(3.20) 


E3,n  VrV^c^VPo^^W172 

1/2    lHv  l' 

E3,y  uy,i  '  Vo  +  si^  <Hy,i>V W  '    -T" 

E-,  u     ,   «  u     „ 

3*2  z,l         z,0 

\  \  '  Un,0  +  bo  Wo^ 
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whe  never 

(3.21)  0<     !i^    <    ^     [i--l],     .0<1. 

o  -o-i 


The 


relations  F     and  E-  are  simply  the  results  of  substituting  p  ■  y"  tin  s 


Icf.   (3.18)1   into   (3.11)  and  (3.12).     The  relation  E,       is  an  immediate  con- 

3,n 

sequence  of  the  longitudinal  part  of  the  flux  conservation  shock  condition, 

B_,  and  (3.7).  E-  _  and  E,   follow  from  (3.5)  and  (3.7).  In  deriving  all 

the  E.-relations  we  have  used  the  fact  that  m  is  positive  in  accordance  with 

the  convention  established  in  subsection  2.1.     The  ratio    |H        /H    |   appearing 

in  the  last  term  of  E~       is  determined  from  Epj   sign  (H  )  appearing  before  this 

term  is  known  since  it  has  been  assumed  that  the  direction  of  this  transverse 

magnetic  field  is  specified.     If  it  is  assumed  from  the  outset  that  H.      is 

tr 

parallel  to  the  y-axis  it  is  only  necessary  to  specify  sign  (H  )  in  addition 
to  (s1-80)/e0. 

In  Figures  U  and  5  we  have  plotted  (p1-pQ)/po  and  (H_/Hx)     respectively, 
versus  s^  0  <  s,  <  5  -  Us     (cf.   (3.20),   the  constant  y  is  assumed  to  be  5/3^ 
for  several  values  of  s   .     These  graphs  will  be  useful  in  Section  5  where  a 
numerical  analysis  of  the  resolution  of  the  Initial  discontinuity  is  given. 


3.3     Fast  gas  dynamical  shocks 

If  we   set  ff.      ,    and  iT      Q  equal  to  zero  in  Eqs.    (2.10)   the  system 
reduces  to  the  system  of  shock  relations  of  one -dimensional  gas  dynamics  having 
a  uniform  flow  transverse  to  the  x-axis.     Assuming  that  (8,-s   )/s     and  the  state 
(p,b,s,u)in  front  of  the   shock  is  prescribed,  we  have  from  the  well- 
known  results  of  gas  dynamics    [see  Courant  and  Friedrichs*-  J,  Chapt.  3ji 

P-i-P  i      2  8i_s~  £.-  "^i  i      2        8,-s 

,  £-2)  .  —^ (-i-i),    -±~±  -    g     I*     (-i_i) 

Ko  1+vl  (s  /e  )         o  Lo  n  +(s  /s  )         o 

(3.22)  l     °  l     ° 

F2  Htr,l     '     hr,0     "     ° 
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(1-M.2)srt      -.1/2 


V  »n,l  *  "n,0  "  bo<— >    U^))  J 


(3.22) 


F3,tr  utr,l     "    utr,0 


ti2+(Sl/s  J       s^-,1/2 


where 


F,  U     -  u    n  -  b         (  i— 2-  )  -2. 

U  n         n,0         o    |_v       ±_^2  Y   J 

H2    »     (y-i)/(y*D  • 
The  fast  gas  shock  by  definition  satisfies  the  relation 


<3-23>  un,l  -  \  >  bl  -   K  Pl^1/2  ■   <  W^  b 


o 


Now  from  (2.7)  and  F,    and  F, 


(3.2U) 


-ln^U/s  J        -.1/2 


u     _   -  U     =  —  (u     „-U    )   =  b,        ±-£-    srt 

n,l         n       Pl       n,0     n'         1    |_       r{l_^)         °J 


The  pure  gas  shock  is  therefore  fast  only  for  those  values  of  (s-,/s  )  for  which 

r  1+H2(S  /B  )     -> 

that  is,  only  for  those  values  of  s,/s  for  which 

8  -      Y-l        l_80  J 


When  s     >  1.  this  relation  is  valid  for  all  positive  values  of  (s, -s   )/b   .     When 
0—  loo 

s    <  1,  the  expression  on  the  right  hand  side  of  (3.26)  is  the  minimum  value  that 

(s..-s   )/s     may  have  in  a  fast  gas  shock.     This  minimum  value  is  clearly  identical 
loo 

with  the  critical  value  of  (s,-s  )/s     in  the  switch-cn  shock. 
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At  this  critical  value  it  is  easily  verified  that  the  values  of  (p ,-p  )/p  , 

(H     -, /H  )    ,  u     ,,  1      ,    and  U     are  the  same  whether  computed  in  accordance  with 
v  y,l'  x     '     n,l'     tr,l  n 

(3.20)  or  (3.22).  It  follows  that  as  the  excess  pressure  ratio  passes  through 
the  critical  value,  the  state  behind  a  switch-on  shock  having  a  given  state  in 
front  passes  continuously  into  the  state  behind  a  fast  gas  shock  with  the  same 
state  in  front. 
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lx.       Slow  and  fast  simple  waves 

U.l  Preliminary  comments  -  Terminology 

Certain  features  of  fast  and  slow  simple  waves  have  already  been  dis- 
cussed in  Section  2.  In  the  present  section,  slow  and  fast  simple  waves  will  be 
treated  in  greater  detail.  Since  the  resolution  of  the  initial  discontinuity  in- 
volves only  slow  simple  vaves,  our  attention  will  be  directed,  in  the  main,  toward 
ascertaining  the  properties  of  such  waves. 

As  already  mentioned, the  flow  adjacent  to  a  region  of  constant  state  is 
always  a  simple  wave  and  the  transition  from  a  zone  of  constant  state  to  that 
of  a  simple  wave  takes  place  across  a  linear  characteristic  [cf.  subsection  2.33  • 
This  characteristic  is  called  the  head  of  the  simple  wave  if  the  fluid  particles 
enter  across  it,  the  tall  of  the  simple  wave  if  the  fluid  particles  leave  across 
it.  If  the  characteristics  in  the  simple  wave  region  of  the  (x,t)-plane  diverge 
with  increasing  time  the  simple  wave  is  said  to  be  a  rarefaction  wave?  otherwise 
the  wave  is  said  to  be  a  compression  wave .  Since  the  resolution  of  the  initial 
discontinuity  involves  centered  simple  waves,  and  since  such  waves  are  manifest- 
ly rarefaction  waves  we  shall  be  concerned  here  exclusively  with  rarefaction 
waves . 

In  subsection  U.2  we  shall  prove  that  in  slow  simple  rarefaction  waves 
the  characteristics  speeds  |u  *  c|  are  at  most  equal  to  the  corresponding  Alfven 
characteristic  speeds  |u  ±  b|,  while  the  opposite  relation  holds  for  these  char- 
acteristic speeds  in  a  fast  rarefaction  wave.  Furthermore,  we  shall  show  that 
centered  rarefaction  waves  are,  in  general,  completely  determined  by  prescribing 
the  state  at  the  head,  and  the  pressure  ratio  at  the  tail.  At  the  head  of  a  slow 
centered  simple  wave  the  transverse  magnetic  field  may  vanish.  Under  these  cir- 
cumstances it  is  necessary  to  prescribe,  in  addition  to  the  state  at  the  head  and 
the  pressure  ratio  at  the  tail,  the  direction  of  the  transverse  magnetic  field $ 
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we  have  already  established  in  subsection  2.3  the  fact  that  the  direction  of  the 
transverse  magnetic  field  is  constant  across  slow  (and  fast)  simple  waves.   If, 
as  we  may  suppose,  the  transverse  magnetic  field  is  parallel  to  the  y-axis,  it 
is  only  necessary  to  prescribe  the  sign  of  this  quantity.  Assuming  the  z-compon- 
ent  of  the  particle  velocity  vanishes  at  the  head  of  the  rarefaction  wave,  we 
shall  show  that  the  transverse  velocity  u  decreases  or  increases  as  we  proceed 
across  the  simple  wave  from  head  to  tail  according  as  the  sign  of  the  transverse 
magnetic  field  is  positive  or  negative.  This  fact  plays  a  role  in  determining 
the  sign  of  H  in  the  rarefaction  waves  which  serve  as  building  blocks  in  the 

y 

resolution  of  the  initial  discontinuity.  In  order  later  to  be  in  a  position 
to  satisfy  the  condition  that  u  vanish  at  the  tails  of  the  rarefaction  waves 
(cf.  (2.21)],  we  shall  find  it  necessary  to  choose  H  to  be  positive. 

The  slow  and  fast  simple  wave  solutions  discussed  in  this  section  were 
first  given  by  K.  0.  Friedrichs  in  [7],  (195U).  In  this  work,  as  already  men- 
tioned in  Section  1,  he  reduced  the  problem  of  obtaining  slow  and  fast  simple 
wave  solutions  in  media  governed  by  the  adiabatic  equation  of  state  to  essential- 
ly that  of  solving  a  single  linear  ordinary  differential  equation  [see  (U.8H. 
It  will  be  shown  that  closed  form  solutions  of  this  equation  exist  for  values 
of  the  adiabatic  exponent,  yt   f°r  which  ©  -  y/(2-y)  is  a  positive  integer.  In 
particular,  a  closed  form  solution  may  be  obtained  when  y  ■  5/3  —  i.e.,  when 
the  equation  of  state  for  the  medium  is  that  of  an  ideal  monotonic  gas.  The  ex- 
pressions for  the  longitudinal  and  transverse  velocity  of  a  slow  (or  fast)  simple 
wave  involve  integrals  which  cannot  be  evaluated  in  closed  form  even  for  the 

_ — — _ 

In  a  recent  work1  J  Owens  has     obtained,  it  appears,   slow  and  fast  simple  wave 

solutions  by  employing  a  procedure  which  differs  somewhat  from  that  of  Fried- 
richs.    His  procedure  is  hased  on  integrating  a  single  non-linear  ordinary 
differential  equation. 
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special  values  of  the  parameter  y  noted  above.  Subsection  U.3  is  devoted  to 
developing  suitable  approximate  expressions  for  these  quantities. 

The  numerical  results  of  the  present  section  are  summarized  graphically 
in  Figures  6  through  12.  Here,  as  in  the  previous  and  succeeding  sections,  all 
computations  are  based  on  the  choice  y  ■  5/3  for  the  adiabatic  exponent. 


U.2  Slow  and  fast  simple  rarefaction  waves 

Slow  and  fast  simple  waves,  by  definition,  correspond  to  the  smaller 
and  larger  roots,  c  and  c„,  respectively,  of  the  equation 

S         X 


(U.l)       ^fftr)2  "  C"2(°2  ~  a2)(p°2  "  ^); 


(U.l)  is  the  result  of  setting  the  last  factor  of  (2.1U)  equal  to  zero.  For 

ryi 

our  present  purposes  it  is  convenient,  following  FriedrichsL  J,  to  introduce 
the  dimensioniess  quantities 

q  -  (c/a) 
(U.2) 

s  -  (pa2)  /  (nH*)  . 

The  quantity  q  is  the  square  ratio  of  the  disturbance  and  sound  speeds  c  and  a 
respectively j  s  is  the  square  ratio  of  the  sound  speed  and  the  Alfven  speed, 

b  -  (^  p-V/2  . 

Since  we  are  dealing  with  a  polytrcpic  gas  undergoing  an  adiabatic  process, the 

2      -1 

square  of  the  sound  speed  a  is  ypP  •  Hence  we  have 

(U.3)       s  -  (yp  /  nH*). 
In  terms  of  s  and  q  (U.l)  becomes 
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(U.U)  (fftr)2     -     (q-l)(s-q'1)H^  . 


Because  the  left-hand  side  of  this  equation  is  non-neeative,  the  first  two 
factors  on  the  right-hand  side  must  have  the  same  sign.     It  follows  immediate" 
ly  that 


(U.5) 


q    <     1,  s  <     —  in  a  slow  simple  wave 


q    >     1,  s>-  In  a  fast  simple  wave 


unless   (H.    ;  vanishes  everywhere  in  the  simple  wave  region,  in  which  case 


(U.6)  q     -     1. 

Under  these  circumstances,  as  may  be  verified  from  an  inspection  of  (2.13)   (the 

D  relations),  the  hydromagnetic  simple  wave  reduces  to  a  one-dimensional  gas 

dynamical  simple  wave  in  which  the  transverse  velocity  u.      is  constant.     Assuming 

that  (ff.    )     does  not  vanish  identically  we  have  from  (U.2)   and  (U.5) 
tr 


(U.7) 


c      <    b    <     a  in  a  slow  simple  wave 

s   *™ 


a  <  b  <  c„  in  a  fast  simple  wave. 


It  follows  immediately  from  these  relations  that  the  characteristic  speeds 
|u  ±  c  |  in  a  slow  simple  wave  are  less  than  or  equal  to  the  corresponding  Alfven 
characteristic  speeds  |u  ±  b | ,  while  in  a  fast  simple  wave  the  characteristic 
speeds  |u  *  cf |  are  greater  than  or  equal  to  the  corresponding  Alfven  charac- 
teristic speeds. 

From  (U.U)  and  the  relations  D0  _  and  D,  of  (2.13)  we  obtain  the  following 
linear  differential  equation  (cf.  [7]  p.  UOj  for  s  as  a  function  of  q: 
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(U.8) 


ds_  9s_ 

dq     +     1-q 


"3 » 

q2(l-q) 


where 


(U.9) 


9     -     Y/(2-Y) 


1  <  y  <  2. 


Now  from  (U.5)  it  is  clear  that  in  a  slow  simple  wave  q  and  s  are  confined  to 
the  region  in  the  first  quadrant  of  the  (q,s)-plane  bounded  on  the  right  by 
q  -  1  and  above  by  the  arc  s-q~,0<q<l    [see  Fig.  6],     Let  (cL...*8*.-)  be 
a  point  on  the  hyperbolic  arc,  i.e.,  let 


(U.10) 


0  < 


<  1. 


tr 


The 


n,  integrating  from  q  to  q^  we  have  from  (U.8) 


©  f 

(lull)       8  -  s(q,q.  :9)  -  -i-  (ri^_)  -  9(l-q)G  J  of  2(l-o>) 


*-(9+l) 


da). 


We  do  not  give  the  corresponding  expression  for  the  fast  waves  since  it  will 
not  be  required  in  subsequent  sections. 

It  will  be  assumed  here  and  hereafter  that  9  is  a  non-negative  integer. 
Several  values  of  the  adiabatic  exponent  y»  1  <  Y  <  2,  for  which  9  is  integral 
are  given  in  the  following  table: 


(U.12) 


Y 

U/3 

3/2 

8/5 

5/3 

9 

,  , ,, 

2 

3 

k 

5 

*  —1  2   -1  1/2 

When  q.      ■  s     ,  then  c  ■  b  »  (nfllP     )       J  that  is,  c  equals  the  disturbance 


speed  of  transverse  simple  waves    [see  subsection  2.5] • 


-  ilO  - 

Note  that  when  the  equation  of  state  for  the  medium  is  that  of  a  polytropic 
monatomic  gas,  then  y  ■  5/3  and  ©  »  5. 

Expanding  the  integrand  in  the  last  term  of  (U.ll)  in  partial  fractions 
we  get 

co  (1-to)  co        k-1     (1-co) 

After  inserting  this  expression  for  the  integrand  into  (U.Li)  we  get,  after 
performing  the  indicated  integration, 

\6      „,.      Af7*\*.      /%r^'^y 


tm-^^<^y^-i 


r 


(U.13) 


-1 

8, 


%avlq    1    «tr  "to     ' 

2     2 
where  q        is  that  value  of  the  square  ratio,  c  /a  ,  of  the  disturbance  speed 

to  the  sound  speed  at  which  cavitation  occurs}  i.e.,  at  which  s  ■  0. 

Let  p_,  b_,  b_  >  0,  H     -i»u     t   and  u     ,  describe  the  physical  state  of  the 
111  y,l       x,i  y,i 

medium  at  a  given  point  of  the  slow  simple  wave.     Solving  for  q  in  (U.U)  and 

2         _2  * 

setting  s  »  s   ,  H.      ■  IT  -,  we  have 

1+E1  +   <*j,l/*x)2  "/l  +  Sl  +   (IV,1/Hx)2  ]"  US1 


(U.lU)  <u  - 

28i 

Since  (q,,g.)  is  necessarily  a  point  on  one  of  the  curves  of  the  family 
8  "  s(q,q.    j©),  where  the  function  s(q,q.rJ©)  is  given  explicitly  in  (U.13), 
we  arrive  at  the  rather  complicated  equation  s,   »  s(q,,q.    ;©)  relating  q.     to 
the  quantities  q     and  s   .     This  equation  may  be   solved  for  q.    ,   graphically  or 


As  already  mentioned,  the  direction  of  H.     is  constant  across  slow  and  fast  simple 
waves.     In  the 
of  generality. 


waves.     In  the  discussion  which  follows  we  may  therefore  set  H    *  0  without  loss 

2 
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employing  appropriate  numerical  methods.     Thus  the  knowledge  of    |H     ,  |   and  s, 
of  a  point  of  a  slow  simple  wave  determines  q.      and  consequently  which  of  the 
8  ■  s(q,q.    ,©)  -  relations  governs  the  dependence  of  s  and  q  throughout  the 
simple  wave. 

Now,  employing  the  equation  of  state  (p/p, )  ■   (p/p.)   ,  we  obtain  the 
following  expressions  for  the  density^ sound  speed^and  disturbance  speed: 

p  •  Pl  (s/Sl)i/Y 
(U.15)  a  -  h(s/Bl)^l)/2r  m  sl/2  bi(8/8i)(Y-D/2Y 

c  -  aq1'2  -  (V-tr)1/2  V^r)1/2  ^l^"1^  • 
From  these  relations  together  with  F0     ,  F„     ,  F0        (cf.   (2.13)]   and  (u.U), 

'»x     '»y     '»* 

we  have 


du       ■     + 
x  y 


1/2 
(U.16)  duy     -     sign  (Hy>1)    [^5]         dux, 


u 


z 


Vl- 


* 


The  factor     sign  (H    ,)  arises  from  the  necessity  of  extracting  the  square  root 

2     -1 

of  the  expression  (H  ,/H  )  ■  (q  -l)(l-qs).  If  we  now  regard  q  as  the  function 

of  s  obtained  by  solving  for  q  in  (U.13)>  q  -  q(e,s  j&),  say,  and  if  we  intro- 
duce the  variable  of  integration 

(U.17)      6"  -  8/str  , 

we  have  on  integrating  from  s/s.  to  s.,/8.  ,  0  <  s  <  s,  <  s. 

*It  will  be  recalled  that  the  direction  of  H  is  constant  across  slow  and  fast 
simple  waves.  Sign  (H  )  is  therefora  constant  across  such  waves* 
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biA 


(U.18) 


where 


x,l         v      s,_    '    IV-V    '  1*  tr'   ' 


0  <  s  <  b1  <  str, 


1 


V  °M  ?  rig"V  ^  ElH«|l     X(s,Sl,8tr),     0<8<Si<str, 


^8>V8tr>  -  <if  ^  / 


Y-l      71' "tr  1        Y+l 


'8/8tr 


5" 

(q/Qt,.)    o-    T   d<r,        o  <  s  <  s1  <  str 


(U.19) 


«*>*i>hJ  '  ^W 


Now,  as  is  easily  verified, 


Y-l       8l/str 


8/6 


•(l-q)/(l-qeav) 

i  -  U/q^  °~ 


tr 


o"      '  d<ri  0  <  6  <  s,  <  s.      . 


(U.20) 


Y-l 

Y-l 

X(8,SliStr)    -    (Ji)    Y  [  X(s/8tr,8tr)    -       ^(-1/«trJ8tr)]     , 

1  '—  J 


where  by  definition 


K(s/s.    ;S_)     :      K(8,s.pJ8tr) 


tr'  tr 


(U.21) 


X(8/str5Str)     '       ^■*8trjetr)' 


We  conclude,  therefore,  that  it  is  only  necessary  to  determine  the  dependence 

of       K(s/s.    Jst    )  and  *.(s/s. _js.    )  on  the  variable  s/s.r  and  the  parameter  s^ 

in  order  to  obtain  the  corresponding  information  about      |^(s,s,js.    )  and 

\.(s,8,;s.    ).     Assuming  this  information  to  be  available  we  have  shown  that  the 

specification  of  the  quantities,  p,  b,  H./H  ,  u     and  u     at  one  point  of  the 

y  x   x     y 

simple  wave  fixes  the  dependence  of  these  variables  on  s/s.   at  all  points  of  the 
the  simple  wave,  in  accordance  with  (U.U),  (U.13),  (U.U*)>  (U.l5),  (U.18)  and 
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(U.19).     We  have  thus  far  implicitly  assumed  that  H    /  0  in  specifying  the  state 

at  one  point.     Now  in  a  slow  simple  wave,  H     vanishes  when  the  disturbance  speed 

c  equals  the  Alfven  speed  b  =   (nIT  p"   )        ,  that  is,  when  q.      »  s"   •     Under  these 

circumstances  in  addition  to  specifying  the   state  (l)  and  b/s.     it  is  necessary 

to  specify  the  sign  of  the  magnetic  field* 

Before  discussing  the  solution  quantitatively,  it  is  worthwhile  at  this 

point  to  note  certain  qualitative  features  of  slow  and  fast  simple  wave  motion 

which  may  be  deduced  without  the  aid  of  a  detailed  numerical  analysis.     In  what 

follows  it  is  assumed  that  the  simple  waves  under  discussion  are  rarefaction 

waves.     We  shall  show  that  in  slow  or  fast  simple  waves  the  density  p,  the  sound 

speed  a,  and  the  disturbance  speed  c  vary  in  the  same  sense  as  the  pressure, 

while  the  magnitude  of  the  transverse  magnetic  field   |H  |  varies  in  the  opposite 

sense.     Moreover,  in  forward-facing  slow  or  fast  simple  waves,  the  longitudinal 
velocity  u     and  the  characteristic  velocity  ux  +  c  vary  in  the  same  sense  as 

the  pressure  while  in  backward-facing  slow  or  fast  simple  waves,  ux,  ux  +  c  and 

the  pressure  vary  in  opposite  senses.     Finally,    we  shall  show  that  in  slow  and 

fast  forward-facing  simple  waves  the  transverse  velocity,  u  ,   varies  in  the 

same  sense  as  or  in  the  opposite  sense  to  the  pressure  according  as  sign  (H ) 

is  positive  or  negative,  while  in  backward-facing  slow  or  fast  simple  wave* 

u     varies  in  the  same  sense  as  or  the  opposite  sense  to  the  pressure  according 

y 

as  sign  (H  )  is  negative  or  positive. 

To  prove  the  first  statement  let  us  rewrite   (U.8)  in  the  following  form: 

||    -     9  q-2(l-q)-1(l-q28). 

Because  of  the  relations  (U.5)  the  right-hand  side  of  this  equation  is  positive. 
Thus  in  slow  or  fast  simple  waves  q  -.aries  in  the  same   sense  as  s.     Since  s  is 
proportional  to  the  pressure    [cf.   (U.3)],  it  follows  that  q  varies  in  the  same 
sense  as  p.     From  (UolU)  we  may  then  conclude  that  in  a  slow  or  last  simple  wave, 
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the  density,  the  sound  speed  and  the  disturbance  speed  vary  in  the  same  sense 
as  the  pressure. 

To  determine  how  |H  |  varies  with  p,  we  write  (U.U)  in  the  form 

H?  -   (q'1-!)  (l-sq)fl?  . 

y  *■ 

Because  of  the  relations  (U.3)  the  signs  of  the  factors  (q~  -1)  and  1-sq  are 
the  same.     Moreover,  the  absolute  value  of  these  factors  and  s  vary  in  opposite 
senses.     It  follows  that    |H    |   and  p  vary  in  opposite  senses  across  fast  and 
slow  simple  waves. 

The  statement  concerning  the  longitudinal,  transverse  and  characteristic 
velocities  follows  from  (U.15)  and  (U.l8)  and  the  corresponding  relations  for 
fast  simple  waves. 


U.3     Approximate  expressions  for    W8/8tr>6tlJ   and  *(s/st,r»s+r) 

We  turn  now  to  the  detailed  examination  of  the  simple  wave  solution. 
In  Figure  6  we  have  plotted     s  »  s(q,qt  jO)    [see  (U.13)]]  against  q,  qcfly  <_  q  <  q^ 
for  several  values  of  q^   .     The  dependence  of    |H  /H |   (which,  from  (U.U),  equals 
Y(q -l)(l-qs)    }.on  s/s     ,  0  <  s/s       <  1,   along  each  of  these  curves  is  depicted 
in  Figure  7.     To  obtain  similar  information  about  u     and  u     it  is  clear,  in 
virtue  of  relation  (U.20),  that  it  is  only  necessary  to  determine  the  dependence 
of     \(s/s     js     )  and  x(B/8trJ8tr)  on  the  variable   (s/str),  0  <  s/8tr  <  1»   for 
the  values  of  s       (-  q~  )  which  correspond  to  the  values  of  q.     employed  above. 
As  mentioned  earlier,  the  complicated  nature  of  the  relation  between  s  and  q 
makes  it  necessary  to  evaluate  the  integral  expressions  for     M(8/Str»5tr^  and 
^•(s/8+   $8+   )  °y  means  of  an  approximation  procedure. 

*  In  all  graphs  appearing  in  this  and  in  the  remaining  sections  0  is  taken  to  be 
5  (see  (U.9)).     The  equation  of  state  of  the  medium  is  therefore  that  of  a 
polytropic,  mon atomic  gas. 


-  U5  - 


Let  us  represent  the  equation  for  the  portion  of  the  8(q,qtr;8)-curve 
through  the  points  (q       ,0)  and  (<ltr>str)  in  the  following  manner: 


U.22) 


(s/str)    |l  *  e(S,str)]   -     (q-qcav),  qcav  <_  Q  <_  ^ 


The  quantity  e(s,s     )  measures,  for  each  q  in  the  specified  range,  the  devia~ 
tion  of  the  s(q,q^   ;©) -curve  from  the  straight  line   segment 

(U.23)  (s/str)     -      (q-qcav)  /   (q^-Q^),  %&v  t  *  <  \T  ' 


Now  the  s(q,q.    ;P)-curves    (see  Figure  6}  are  almost  linear;  we  may  therefore 
neglect  e(s,£.    )  without  incurring  too  great  a  loss  in  accuracy.     In  the  analysis 
which  follows  (U.23)  will  be  employed  in  lieu  of  (U.22). 
It  will  now  be   shown  that 


He/stristr)  -  5(qcav/qtr) 


1    - 


*  -(•/■*>*    >+fr   {l-(*/°tr> 


(U.2U) 


11 


W  ^    •  (S/Str»: 


,     0     <     S/S  <     1, 


,(    ,  x        f^A     1       a   \       6+8  _,6     1.     s   ,    fc  362+268-p2  D,ll     1       s    , 

V(s/3tr,str)   -    |B(F,  j,  _)  -  -gE.  B(5,  ?;  —)  +         -     PP     B(T,  ?,   — ) 


(U.25) 


tr 


B  6-8^-386  -56^  D/16     V,   .     ».        1     o  «. 
+  E £-^E B(T,   2-j   str)+   ...J,  0  < 


s/8tr  -  lu 


Here 


(U.26) 


B(p,q;s/str) 


c-p-^i-o-)^1  d<T, 


s/s 


tr 


and  the  quantities,  a,  8  and  6  are  defined  as  follows: 

*  The  quantities  p  and  q  employed  here  are  in  no  way  related  to  the  pressure 

2     2 

and  the  square  ratio  q  «  c  /a     . 
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°     =      (qtr-qcav)   I  qcav 
(U.27)  P     «     Kr-%^  I  ^"W 

5    "     (qtr-qcavJ  '  \r     • 

To  obtain  the  above  expression  for      |^y(s/s.    ;s,    )  from  the  integral  ex- 

1/2 

pression  for  this  quantity  [cf.  (U.19)3  **  expand  (q)  '     in  powers  of 

(q-q   )/q   .  substitute  for  q-q    from  (U.23);  then  we  multiply  through 
cav   cav  cav  '  r 

by  O"  "   '     and  integrate  from  s/s   to  unity.  To  obtain  the  corresponding  ex- 

tr 

pression  for  X(s/s     ;s     )    (cf.   (U.19)]]  we  first  substitute  for  q  in  (l-qs)"  ' 
from  (U.23)  to  get 

1/2    |-  -,-1/2 


\l  -<8/8tr>]"  \_  ~  5<s/str>] 


1/2 
and  expand  the  last  factor  in  powers  of  6(s/s      ).     Then  we  expand  (i-q)         in 

powers  of  (q-q       )/(l-q       )  and  express  q-q^-  in  terms  of  s/s.      by  means  of 

(U.23).     The  resulting  power  series  is  then  multiplied  by 

r-  1/2 


•  -U/5 


[l  -  (s/str)J 


r  i_1/2 

and  by  the  power  series  expression  for     1  -  p(s/s.    )  and  then  integrated 


termwise  from  (s/s.    )  to  unity, 
tr 


In  Figures  8  and  9  we  have  plotted  a,   6  and  p  versus  q.    ,  0  <  q.     <  1. 
The  limiting  values  of  these  ratios  near  q.      -  0  are  obtained  from  the  fact 


that* 


and 

1  U 

^^  ^cav'-^^-^tr^     <kr*1- 

*~ 

bee  Appendix  I. 
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It  should  be  observed  that  a  and  6  are  less  than   .2  in  the  range  0  <  q.     <  1 
T-rhile  (3  is  less  than  .5  for  0  <  q.      <   .975  but  increases  rapidly  to  unity  when 
.975  <  q.     <  1«     The  series  expression  in  the  right-hand  side  of  (U.23)  there* 
fore  converges  rapidly  for  all  q.    ,  0  <  q.     <  1;   the   same  is  true  of  the  right- 
hand  side  of  (U.2U)  over  the  smaller  q^   -interval  0  <  q.      <   .975. 

It  is  desirable  to  rewrite  the  expression  for  \(s/s.    ;s,  J  somewhat,  in 
order  to  put  in  evidence  the  relative  contributions  of  the   individual  terms  of 
the  sum  and  to  have  an  expression  better  suited  for  computational  purposes.     To 
this  end  let  us  write  B(p,qjx)  in  the  form 


(U.30)  B(p,qjx)     -     B(p,q-,0)    L  -  Ix(p,q)l 


}  0    <    X   <    1      j 


where 


(U.31)  I  (p,q)     " 


Note  that 


) 


x 


^^""(i-o-J^do- 


B(p,qjO) 


(U.32)  0  <  Ix(p,q)  <  1,  0  <  x  <  1   ^ 

and  that 

1 

(u.33)         B(p,qjo)  -    I  o-p-^i-o-^-V  -  HzlEM, 

I  T(p-q) 

the  last  equation  being  a  statement  of  the  well-known  relation  between  the  gamma 
function  and  the  beta  function.     Substituting  for  B(p,q;x)  in  (U.2U)  and  making 
use  of  the  fact  that 

B(6/5,    1/2;  0)     -     (2/7)  B(l/5,  l/2j  0) 
(U.32)  B(ll/5,  1/2j0)     -     (2A)(i2/7)  B(l/5,   l/2;0) 

B(l6/5,  l/2;0)     -     (2/7)(l2/17)(22/27)B(l/5,   l/2;0) 
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we  have  finally 


\(s/s.    js,    ) 
v   '   tr     tr 


B(l/5,l/2jO) 


(U.35) 


>2  03 


;3v2v/l2w22. 


15 


T 


-,'? .  5> 


+    ... 


where 

c  -,         rd/5)r(i/2) 

(U.36)  B(S,  UO)     '     -  6.286. 

5  2  r(7/io) 

When  p  >  6/5,  numerical  values  of  I    (p,l/2)  for  values  of  x  in  the   interval 
0  <  x  <  1  may  be  obtained  from  K.  Pearson's  tables  of  the  incomplete  beta 
functions L <J by  interpolating  between  the  tabulated  values  of  p.     The  first 
term 


B(l/5,l/2jx)     -     B(p,qjO)fl  -  Ix(p,q)l 


-1/2 
may  be  obtained  directly,  for  example,  by  expanding  the  terms  (l-o~  )     '     as  a 

power  series  in  o~    and  integrating  termwise. 


The  dependence  of     K,(s/8tr»stI.)  and  X(s/8tr*Str^  0n  s^str»  °  -  s/\r  -  1» 
for  q^     -  s"1  «  0.05  and  q.      -  0.975  is  depicted  in  Figures  10  and  11.     The 
curves  corresponding  to  q^     lying  between  these  values  are  not  shown.     Figure 
12  is  added  in  order  to  facilitate  interpolation  between  these  curves}  in  it 
we  have  plotted    r\(0js     )  and  \(0jstr)  versus  q^,  0  <  qtr  <  1.     This  informa- 
tion is  useful  for  interpolation  since  it  is  found  that 
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K[s/Str>str]  "  ^/^.OS)-1]    K|o;«tr]  -  ^(.OS)"1] 

— __ ____  ^ _ 

(U.37) 

X[S/StrjStr]  -  ^/^tr^^)-1]     *|o,stp]  -  ^(.OS)"1] 


to  a  sufficient  degree  of  accuracy. 

From  the  integral  expressions  for   '^s/strjstr^  and  ^8^etrJStr^  it  iS 
easily  verified  that 


^s/8tr'str)=l^(0}Str)     "     5  V^f  ^stv)VS>     W*tr>  *  ° 

(U.38) 

\.(s/strjstr) «  X(Osstr)     -     5  U/str)1/5,  (s/str)->0. 

These  asymptotic  results  are  independent  of  the  linear  approximation  employed 
in  deriving  the  series  expression  for      K(s/str}str)  and  \(s/str;str). 

To  prepare  for  the  analysis  of  the  next   section  it  will  be  useful  to 
obtain  the  corresponding  asymptotic  relations  for  ^s/str'fetT^   and      *\WstrJStr^ 
as  s/s       approaches  unity.     In  the  neighborhood  of  s/str  -  1  the  sCqjq^jG)- 
curve  may  be  approximated  by  the  straight  line, 

^    '  (8-8+ J. 


q     "     ^  +  di  1  (6"Str 

l^^r 


Employing  (U.31)  we  have,    since  q^  ■  a^r   , 


V  -  ¥  [l  -  C/'tr'] 


It  follows  that 


i-q.  =r  *£.    (i  -  (8/str)]  Ji  -  sk  £  -  (»/»tP>] 
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Thus,  retaining  only  first  order  terms,  we  get 

(U.39)        1  -  qs  S  §gi  [l  -  (s/str)J  . 

Similarly,   neglecting  higher-order  terms,  we  have 
1   -  q     S     i  -  qtr. 

Employing  these  approximations  in  the  expression  for  X(s/s     js     )  we  have 

tr     tr 


.(s/s+    ;s+    )   -    (l-q        )"1/2    f  [i=^f  a7  %T  do- 

tr*   tr'  ^cav'  J  |_I-q8  _ 


A 

s/str 

1 


^cav  ,  . 

S/B. 

'   tr 


Integrating  and  retaining  only  higher-order  terms,  we   get  finally,   since 
6/6+1  -  y/2    (cf .   (U.9)], 


1  1 

(U.U0)         X(s/str;str)  a    jiY  (l=^  )J?  ji  -  (s/str)]r,         s/s^  -  1. 

The  corresponding  expression  for     K (s/st_;s.    ), 

(U.U1)  H*l\r\\J     '     T1  '  (S/Str}]  ' 

follows  immediately  from  the  definition  of  ^\(s/ptr»str)» 
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5 .       p.3 solution  of  the  shear  flow  discontinuity 

5.1       Preliminary  considerations;  notation 

As  indicated  in  Section  1,   the  resolution  of  the  initial  discontinuity 
may  be  effected  by  two  fast  shocks  of  equal  strength  moving  away  from  the  initial 
discontinuity  with  equal  but  opposite   speeds,   followed  by  two  slow  centered  rare- 
faction waves.     These  slow  centered  rarefaction  waves  are,   like  the  shocks, 
symmetrically  disposed  with  respect  to  the  origin  at  each  instant,   and  connect 
the  constant  states  behind  the  shocks  to  an  expanding  zone  of  constant  state  be- 
tween the   simple  waves  or  to  a  vacuum  zone  which  expands  with  an  appropriate 
"  escape   speed"     according  as  cavitation  does  or  does  not  occur „     It  will  be  re- 
called that  the  fast  shocks  could  either  be  fast     pure  gas  dynamical  shocks  or 
switch-on  shocks,  and  that  the  choice  of  these  waves  was  dictated  by  the  absence 
of  a  transverse  magnetic  field  in  front  of  the  shocks.     In  light  of  the  analysis 
of  Section  3  we  now  know  that  the  resolution  of  the  initial  discontinuity  is 
effected  by  fast  gas  shocks  unless  s   ,  the   square  ratio  of  the  sound  speed  and 
the  Alfven  speed  in  the  undisturbed  medium, is  less  than  unity  and  the  shock 
strength,  measured  in  terms  of  the  pressure  ratio,  is  below  a  certain  critical 
value,  in  which  case  the  resolution  is  effected  by  switch-on  shocks. 

In  further  specifying  the  nature  of  the  resolution  it  is  assumed    [see  sub- 
section 2*14]  that  (i)  the  transverse  component  of  the  particle  velocity  is  parallel 
to  the  y-axis  and  that   (ii)   the  pressure,   density,  and  transverse  magnetic  field 
components  are  even  functions  of  x,  while  the  velocity  components  are  odd  functions 
of  x.     Assumption   (i)  implies  that  the  transverse  component  of  the  magnetic  field 
is  parallel  to  the  y-axis    [see  subsection  2, IT),     Assumption   (ii)  implies  the 
vanishing  of  the  transverse  and  longitudinal  velocities  in  the  zone  of  constant 
state  behind  the  rarefaction  waves,  assuming,  of  course,  that  cavitation  has  not 
been  produced.     In  case  the  magnitude  of  the  initial  shear  flow  discontinuity  is 
large  enough  to  produce  cavitation, the  vanishing  of  the  velocity  behind  the  rare- 
faction waves  must  be  replaced  by  the  condition  that  the  electric  field  vanishes 


-  52  - 

at  the  tails  of  the  rarefaction  waves    [see  subsection  2.U,  equations  (2.22)- 
(2.30)T)  •    Under  these  circumstances  the  tails  of  the  rarefaction  waves  are 

allowed  to  move  in  the  y-direction  with  equal  and  opposite  transverse   sreeds 

it  ■ 

and  in  the  longitudinal  direction  with  equal  and  opposite  escape  speeds  ♦ 

The  object  of  this  section  is  to  prove  that  this  motion  is  indeed  a 
solution  of  the  initial -value  problem  formulated  in  Section  1  and  to  provide 
the  basis  for  a  complete  description  of  the  motion.     Before  doing  this,   however, 
it  will  be  convenient  to  establish  here  the  notation  to  be  employed  and  to 
describe  briefly  the  procedure  to  be  followed  in  subsequent  subsections. 

Because  of  the  symmetry,  it  is  clear  we  need  only  consider  the  motion  in 
the  region  x  >  0.     In  Figure  3a  we  have  depicted  the  basic  flow  regions  in  the 
positive  half  of  the   (x,t)-plane.     These  basic  flew  regions  are: 

Region   (0)    :  The  region  of  constant  state  in  front  of  the  forward-facing 
shock; 

Region   (l)   :   The  region  of  constant  state  behind  the  shock; 

Region  (2)   j   The  flow  region  of  the  slow,   centered,  simple  forward-facing 
rarefaction  wave} 

Region  (3)    s  The  region  of  constant  state  following  the  tail  of  the  forward- 
facing  simple  wave. 

Similar  quantities  in  the  various  regions  will  be  distinguished  by  appending 
appropriate  subscripts.     The  state  at  a  point  in  region  j,   j   ■  0,   1,   2,   3>  will 
be  characterized  by  the  values  p.,  p.,  H      .,  H  ,  u      .,  u      .  or  alternatively 
by  the  values 


s 


(-  YP3*4),  bjC-Vf^pJ1  ),   (Hy>j/Hx),   (u^j/bj)  and  (u^/ty  or  C^j/.j). 


It  will  be  recalled,  that  p     is  the  pressure,  p.  the  density,   R^  the  fixed  longi- 
tudinal magnetic  field  component,  H     .,  the  transverse  magnetic  field  component, 

b1  the  Alfve'n  speed,   a.  the  sound  speed,in  Region  j,   J  ■  0,  1,2,3  and  s.(-(a./b.)   ) 
"  «  3  3     3 

the   square  ratio  of  the  sound  speed  to  the  Alfve'n  speed  in  Bfcgion  j,   j  -  0,1,2,3. 
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The  complete  description  of  the  motion  entails  the  knowledge  of  how  the 

physical  state  in  regions  (0)-(3),  the  velocity  of  the  shock  front,  and  the 
head  and  the  tail  of  the  simple  wave  depend  upon  the  thermodynamic  state  in 
front  of  the  shock,  on  H  ,  and  on  the  magnitude  of  the  initial  discontinuity 
in  velocity  2u  0.  To  establish  this  dependence  we  recall  that  the  state  be- 
hind a  fast  shock  and  the  state  at  the  tail  of  the  rarefaction  wave  depends 
only  on  the  pressure  ratios  P-./P  (■  s  /s  )  and  P^/P-,  (*  so/s-, )  respectively, 
assuming  the  state  in  front  is  specified  [see  Sections  3  and  U]  .  Since  the 
state  (l)   behind  the  shock  is  the  state  in  front  of  the  rarefaction  wave  it 
follows  that  state  (j)  at  the  tail  of  the  rarefaction  wave  depends  upon  the 
two  parameters  p./p  and  P-i/Pi  in  addition  to  the  state  in  front  of  the  shock, 
If  cavitation  has  not  been  reached,  then  p_  •  s  f   0 .  From  the  expressions 

for  u  and  u   [cf .  (U.l8)-(U.19 f)  in  the  simple  wave  region  we  obtain,  after 

■*■  y 

setting  u     ■  u     »  0  at  the  tail,  two  relations  of  the  form 

y   x  ' 


Ri  [JW'  (S1/S0)J  ^(A^o]  "  °> 


1,  2. 


By  a  combination  of  the  graphical  and  numerical  methods  we  then  obtain  the 

dependence  of  8.,/s,  and  s,/s  on  u  „/b  for  several  values  of  the  parameter  s  . 
3  1     1  o     y,0'  o  ^         o 

It  can  be  shown  that  cavitation  is  produced  when  the  ratio  u  „/b  exceeds 

y  y,0'   o 

a  certain  minimum  value  depending  on  the  parameter  s   •     In  this  case   s_   («  'YpVuir) 
vanishes  at  the  tail  of  the  rarefaction  wave.     This  condition  in  conjunction  with 
the  second  condition  of  (2.22)  leads  to  a  relation  of  the   form 


sl.W^o/V'   So]     "     ° 


from  which, again  by  numerical  and  graphical  means,  the  dependence  of  the  ratio 
s,/s^  on  the  ratio  u  n/b„  and  on  s„  is  established.  It  should  be  observed  that 

i    o  y>o'  o  o 

the  state  in  front  of  the  advancing  shocks  is  characterized  by  the  values  p   , 

u„  r\P°^i   s   .  b   .     It  is  therefore  fortunate  that  the  relations  R.    ■  0,  i  ■  1,2, 
y,0     o       o       o  i         '  *   ' 

and  S  -  0     do  not  depend  on  p     and  b     explicitly. 


We1 leave  open  for  the  moment  questions  involving  the  direction  of  the  trans- 
verse magnetic  field.   See,  however,    subsections  U.l,   5.21,  and  5.31» 
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In  subsection  5.2  we  summarize  the  formulas  giving  the  state  behind  the 
fast  shocks  in  terms  of  the  pressure  ratio  Pt/P0  ■  si/s0  and  the  state  in  front. 
In  subsection  5.3  the  corresponding  expressions  characterizing  the  states  in  the 
simple  wave  region  are  given  in  terms  of  the  ratio  P9/P-1  ■  so/si  anc*  *"he  one- 
parameter  family  of  states  in  Region  (1)  at  the  head  of  the  rarefaction  wave. 
The  program  for  establishing  the  depandence  of  s,/s  and  s_/s,  on  u  0/bQ  for 

various  values  of  8  is  outlined  in  subsection  5.U.  The  results  are  summarized 

o 

in  the  Figures  l£<:  and  16.  In  expressing  the  dependence  of  s-,/s  on  the  state 

in  the  undisturbed  medium  it  is  also  useful,  for  reasons  which  will  be  apparent 

later,  to  plot  s,/s  versus  u  ~/a  ,  where  a   ■  b  Vs~      is  the  speed  of  sound 
'         *  i'  o        y,0'  o*       o     o  y  o  I        ^ 

in  the  undisturbed  medium.  The  curves  corresponding  to  s  >  1  are  shown  in 

Figure  15a,  those  corresponding  to  s  <  1  in  15b.  The  limiting  forma  of  the 

dependence  of  s,/s  ,  8o/s,  (IL.  o/H  )  and  other  quantities  of  interest  on  u  0/b0 

corresponding  to  small  and  large  values  of  u  n/b  are  discussed  in  subsection 

5.5«  On  the  basis  of  the  results  of  this  and  the  preceding  subsection  it  is 

a  routine  matter  to  plot  the  remaining  state  variables  in  regions  (l)  -  (3) 

versus  u  n/b  or  any  related  quantity.  For  the  sake  of  brevity,  however,  we 

shall  merely  plot  the  ratio  (H  -/H  )  versus  (u  ~/a  )  and  (u  »/b  )  and,  in 
*  r  y,3  x'         y,0'  o'      y,0'  o'    ' 

addition,  (u  _)   /b  and  (u  -,/b  )  versus  u  _/b  .  Here,  (u  ,)    is  the 
'  v  x,3  esc'  o     v  y,3  o'        y,0'  o      '   x,3  esc 

longitudinal  w  escape  speed"  and  u  ,  the  transverse  velocity  at  the  tail  of 

y>  j 

the  rarefaction  wave  when  (u  n/b  )  is  large  enough  to  produce  cavitation. 

y,u     o 


5.2     The  state  behind  the  fast  shocks  assuming  the  pressure  ratio  and  state 
in  front  is  prescribed 

5.21    The  switch-on  shock.     Let 

pl  sl 

(5.D  1  -    ±    -    ±  . 

*o  o 

From  (3.20)  we  have,  on  taking  into  account  the  conditions  in  front  of  the  switch- 
on  shock  and  the  fact  that  the  shock  is  forward-facing, 
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s. 


prpo 


Y-l 


.1  ^  i  +  illzi)  (  ^i, 


1/2 


rHy,i)2 


*<¥*> 


l,Y-l   /l"Po^ 
( V  ^~p^~)  +   So 


(5.2) 


'3,x 


'3,y 


u 


x,l 


u. 


y,i 


u 


y,o 


yA 


u 


o  v  pl 


whenever 


(5.3) 


v-.  v%a-o 


s^  <  1, 
o 


Note  that 


(5.U) 


ux,l  +  bl 


Ux» 


that  is,  the  characteristic  Alfven  speed  behind  the  shock  is  equal  to  the  shock 
speed 

It  should  be  observed  that  the  negative  sign  has  been  employed  in  front  of 
the  term  involving  the  transverse  magnetic  field   (cf.  E^     ,   (3.20)].     In  making 
this  choice  of  sign  we  are  anticipating  the  fact  that  the  switch-on  shock  is 
followed  by  a  rarefaction  wave  of  the  type  discussed  later  in  subsection  5.31« 
It  will  be  seen  that  H     is  positive  in  this  rarefaction  wave.     In  order  to  make 

y 

a  continuous  connection  at  the  head  of  the  rarefaction  wave  with  the  state  in 
front,  H     ,  must  also  be  chosen  positive.     Equation  S_       is  a  consequence  of 
this  requirement  and  the  fact  that  the  shock  is  forward-facing. 


-  56  - 

5.22  The  fast  gas  shock.  From  (3.22)  we  have,  on  taking  into  account 
the  conditions  in  front  of  the  shock  and  the  fact  that  the  fast  shock  is  forward- 
facing, 


G 


'l-o     _  1-jT 


i         T-  -rx  K  <  ""  —, 


(v?  -«J     -V1    "     T*-    (V7-L)'     r"p"1 


o 


i+H-  >7  o  u-  +v/ 


°2  Hy,i    -    Hy,o    -    u 


2  X 


(5.5)  LY(t4+>()J 


G3,y  uy,l    "    uy,0 


2  1 

.(nS>7)8 ol  7 

G.  U       -     b 

U  x  o 


L  YdV)  J 


whenever 


(5.6)  so    >    1 


or  whenever 

(5.7)  8Q    <     1    but    >j-l    >     |Xj.     (s^1  -1). 

Here 

(5.9)  u2    -    &    .       § 

From  (3.23)  we  have 

(5.9)  uxA     +     b±    <    Ux  ; 

i*e.,  the  Alfve'n  characteristic  velocity  behind  the  shock  is  less  than  the 
velocity  of  the  shock  front. 
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In  the  resolution  of  the  initial  discontinuity,  gas  shocks  will  be  employed 

when  s     >  1.     When  s    <  1,   switch-on  shocks  or  gas  shocks  will  be  employed 
o  —  o 

according  as    ri  -1  is  less  than  or  greater  than  the  critical  value     2y/(y-1)   (SI  -1)> 
As  shown  in  Section  3  the  state  behind  the  switch-on  shock  is  identical  with  that 
behind  the  pure  gas  shock  when     V^ -1  -     (2y)/(y-1)   (*£  -1)»   s0  tnat  as     ^ -1  in- 
creases through  the  critical  value  the  switch-on  shock  changes  continuously  into 
a  fast  pure  gas  shock. 

5.3  The  state  in  the  simple  wave  region  (Region  ,{2))   behind  the  shock 

5.31    The  RS  combination*   fsQ  <  i;    yj  -i  <  -j  tZ-  (s^  -1)  M    .     Since 
the  state  in  front  of  the  rarefaction  wave  is  the  state  bshind  the  switch-on 
shock  and  since  the  rarefaction  wave  is  forward-facing  in  the  .region  x  >  0  we 


have  from  (U.15),   (U.18)  and  (5.2) 

RS 


1  P0  Po        Sl  -    Z     1  - 


RS2  (Hy,2/Hx^     "     (qi1-D(i-V2)'         Hy,2>° 

11  11 

RS3,x  Tf  '   |Vpo>?  "  <P<Al>?]  ~  *  ^1^  Vst/  \<82>  Vstr> 

K  ,    .      (p«/pJ     r  -r? 


RS3,y 


o  o  x  —  -1 

1  1  Y-l  1  1  1  Y-l 

8i\7q2  \2/s2v    k  /'po\5/si\^q2\  Yjy^ 


where  exact  and  approximate  expressions  for    K^Sg^js^)  and  Xts^s-^s^)  are 
given  in  Section  k   (cf.   (U.19),   (U.20),   (U.2U)  and  (U.35fJ  and  where 


*  The  expression   "  the  RS-combination "    should  be  read  as   "  the  rarefaction  wave- 
switch-on  shock  combination M  • 
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P, -P.  s 

(5.11)  i-2. 

0 


^[.i.|i.^<^|^ 


[cf.  equation  S  ,    (5.2)3. 

At  the  as  yet  undetermined  head  of  the  simple  wave,  s„  ■  s, ,   q?  =  q  ,   and 


'5.12)  (H^/y2  -   <HyyH/  -  2(-i-^)    |ll-s0,  -  t£  (-il£)] 

(cf.   S2,    (5.2)].     Solving  for  qx  in  RS1,   (5.1C),  we  get 

(513)  q     .    <^»l>    ''"7)1/"x'2  '/iK  *    Cy.l^x'^2   -  ""l 

1  2=1 

where    [see   (5.1)3 

(5.1U)  s1     =     Vf  sq   . 

Assuming  for  the  moment  that  s     is  fixed  we  may    regard      (5.13)   and  (5. Hi)   as  the 
parametric  representation  of  a  curve  in  the   (q,s) -plane,     >7  serving  as  the  para- 
meter.    In  Figure  13  we  have  plotted  sx  versus  q,    for  several  values  of  s    .     The 
points  on  these  curves  lie  below  the  hyperbolic  arc  s  *  q"   ,    .2  <  q  <  1.     Only 
at  the  critical  value     >f    ■  1  +  -X-  (s~  -1)  is   (q,,s,  )  a  point  on  this  arc,  and 
by  definition    >1     is  less  than  this  value  in  a  switch-on  shock.     We  are  therefore 
led  to  the  relation 

(5.15)  0     <     Sj/s^.    <     1 

which,  in  virtue  of  RS.  ,  implies  that 


(5.16)  c2  <  bl  J 


i.e.,  the  disturbance   speed  in  the  simple-wave  region  is  less  than  the  Alfven 
disturbance   speed  behind  the  switch-on  shock.     Combining  this  result  with  that 
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stated  in  (5.U)  we  arrive  at  the  fallowing  relation  between  the  velocity,  U  >  of 
the  shock  front  and  the  characteristic  speed,  u  ,  +  c  ,  at  the  head  of  the  rare- 

A  y  X  X 

faction  wave: 

(5'17)  ux,l  +  cl  <  ux,l  +  bl  "  Ux  • 

Thus  the  rarefaction  wave    "  follows"     the  switch-on  shock. 

It  should  be  observed  that  sign  (H     .  )  has  been  chosen  positive  in  passing 
from  the  expression  for  u     in  (U.l8)  to  that  given  in  RS-      .     This  choice  of  sign 
is  made  to  meet  the  requirement  that  u     _  vanish  at  the  tail  of  the  rarefaction 
wave  for  all  values  of  u    .  which  do  not  exceed  the  minimum  value  necessary  to 
produce  cavitation. 


5.32     The  RG  combination*  Is     >  1  or  s    <  1  but  v?  -1  >  ^X-  (s^-l)  1.     At 
the  head  of  the  rarefaction  wave  in  the  RS  combination  v?e  have  q^  ■  q_    and  s_  ■  s,, 
It  will  be  shown  below  that 

(5.18)  qx  -  qtr,  s1  -  str     -     q£   . 

From  (5.5),  which  expresses  the  state  behind  the  fast  shock  in  terms  of  the  ratio 
ri     ■  p./p     and  the  state  in  front,   and  from  (U.15)  and  (U.l8),  which  expresses 
the  state  in  the  simple  wave  region  in  terms  of  the  ratio  6^/8.      and  the  state  at 
the  head,  we  then  have 

Po         r..2+ 


Ko  '—  l+|i  *\  -1       tr 


(5.19) 


RG2  (Ky,2/Hx)2     '     (q21-1)(l-^2S2) 


The  expression  "the  RG  combination"  should  be  read  as  "  the  rarefaction  wave- 


pure  gas  shock  combination"  . 
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1  1 


uy,2 


RCVx        -tr~  s  ^  "^ 


3,x  bo 


2X.    -, 


"(1-H')8U 

_y(^2+^  ) 


x    -....2._    -,? 


r 


[^t ]  HWt,> 


■  H   +VJ 


(5.19  cont'd)  1 


RG 


Etr^cavO*  ri+  2, 


Uy,2       uy,0  ^trv"  Hcav<-1      j   1+g.-^ 

3»y  ~^~~       b      "  7  ~2~~~ 


-T2- 


^W^ 


T±  Y-i 

1/2  2y^  p     1/2  2v 

Mcav  tr  1  nr       tr 

Exact  integral  expressions  for    ^(62/str,str^  and  X^s2//str,Str^»  °  -  s2/str  -  l* 
may  be  obtained  from  (u.19).     Approximate   series  expressions  for  these  functions 
are  plven  in  (U.2U)   and  (U.35).     Plots  of    Y\    and  K  versus  s„/s     ,  based  on  these 
approximate  series  are   given  in  Figures  10  and  11. 

That  q«  ■  q.      and  8     ■  q~     at  the  head  of  the  rarefaction  wave  is  a  conse- 
quence of  the  fact  that  H     ,    vanishes  behind  the  gas  shock    (cf.  equation   G0,    (5.5)1. 

y,j.  •—  c  -» 

At  the  head  of  the  rarefaction  wave  we  then  conclude  from  RG,,  that  q.  «  1  or 

q„  "  s^  .  It  is  easily  verified  that  when  q2  ■  1,  the  hydromagnetic  simple 

wave  reduces  to  a  pure  pas  dynamical  wave.  Such  a  simple  wave  cannot  be  employed 

to  effect  the  resolution  of  the  initial  discontinuity  since  u  „  (»  u  „)  remains 

y,2    y,C' 


constant  across  the  weve .  The  remaining  alternative  is  therefore  q?  ■  q.   and 
s_  ■  s   -  q~  at  the  head  of  the  rarefaction  wave. 

From  RG.  we  conclude  that  c2  "  b,  at  the  head  of  the  simple  wave.  We  have 
therefore  verified  the  fact  that  the  head  of  the  simple  wave  follows  the  fast 
gas  shock  (see  (5.9)]. 

Note  that  we  have  again  chosen  sign  (H  )  to  be  positive  in  passing  from  the 

expression  for  u  in  (U.18)  to  that  given  in  RG,  ,  (5.19).  As  is  the  case  with 

the  RS  combination  this  choice  of  sign  is  necessitated  by  the  requirement  that 

u  _  vanish  at  the  tail  of  the  rarefaction  wave  for  values  of  u  -  below  the 
y>3  y,0 

minimum  value  required  to  produce  captation. 
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5.U  Program  for  establishing  the  dependence  of  ~n(=  s,/s  )  on  (u  n/b  ) 

\     J.  o       y  j  u  o 

for  several  values  of  s  . 

o 

5.U1  Preliminary  remarks.  Our  objective  is  to  obtain  a  solution 
satisfying  the  conditions 

(5.2C)  ux>3  -  0,        uy>3  =  0 

in  the  region  of  constant  state  behind  the  simple  wave  unless  cavitation  occurs 
in  which  case   (5.20)  is  replaced  by 

(5.200  P3-0,  uy>3Hx  -  uXj3Hy>3     -     0 

at  the  tail  of  the  rarefaction  wave.     We  then  investigate   the  dependence  of 

this  solution  on  the  initial  state.     This  initial  state  will  be  characterized 

bv  the   quantities  u     „/b   ,   p   ,  b   ,    s   .     As  an  inspection  of  (5.10)  and   (5.19) 
J  ^  y,C     o       o       oo 

reveals  that  it  suffices  to  know  how  ri    and  3.,/s..   depend  upon   (u     _/b   )   and  on 

the  initial  thermodynamic  state   (p   ,b  ,s   ).     Actually,   as  will  soon  be  apparent, 

Y>       and  s  /s     depend  only  on  (u     Q/b   )   and  s   .     Consequently  the  amount  of 

computation  that  is  necessary  to  establish  the  dependence  of  the   solution  on  the 

thermodynamic   state  in  front  is  materially  reduced. 

The  procedure  for  establishing  the  dependence  of    vj     and  s./Sj^  on  (u    q/^>q) 

and  s     falls  naturally  into  two  parts  according  as  cavitation  has  or  has  not 

occurred.     Let  us  first  prove  that  cavitation  always  occurs  when    v^    is  large 

enough.     For  fixed  p   ,   s   ,  b     let     vj  be  the  minimum  value  of  the  ratio 

*>  "o       o       o  ( cav 

s,/s     (=  p,/p   )  at  which  cavitation  occurs.     For  this  value  of    vi      u         still 
1'    o  vL     o  xiJ 

vanishes  at  the  tail  of  the  simple  wave.     Solving  RG3  x,    (5.19)   for     vjcay,   after 
first  setting  the  left-hand  side  of  this  equation  equal  to   zero,  we  obtain* 


(5.21)  ^    ■    <**>»/<*»»>' »«^>      ,      A.[K<o,v)2/5s0], 


The  quantity    yn    is  at  least  unity;   this  fact  requires  the  use  of  the  positive 
sign  before  the  radical  in  (5.21). 
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for  all   s     such  that 
o 


(5.^2)  V]  -  i     >     £l_     (s"1-!). 

w       '  I cav  —    y~1         ° 

Note  that  K^(0;s  )  versus  ck  »  s~  has  been  plotted  in  Figure  12;  it  should  be 

observed  that  ^(0;s.  )  varies  only  slightly  with  s.  .  When  s  is  fixed,  the 

tr  x»r  o 

right-hand  side  of  (5.21)  therefore   varies  only  slightly  with  s     .     The  quantity 

T(s.    )        /s  1  «     vi  may  then  be  obtained  by  a  rapidly  convergent  iteration 

LV   tr  cav'   oJ  'cav      J 

procedure.     In  Figure  1U  we  have  plotted    >i  versus  s    .     This  curve  is  inter- 

sected by  the  curve     *]  -1  -     (2y)/(y-D   (s~  -1)   at  sq  ^0.05.     When  sq  is  below 
this  value  the  expression  for     r>  piven  in  (5.21)  must  be  replaced  by  a  similar 

expression  which  results  from  SO..     .   (5.10),   after  setting  u     0  »  0  and  solving  for 

j,x  x,£ 

vo    •  Since,  however,  values  of  s  below  .0625  will  not.  be  employed  in  our  com- 
(  cav  o 

putations  it  will  not  be  necessary  to  present  this  result  here.  It  is  neverthe- 
less worthwhile  mentioning  that  in  this  case  too  ^__v  depends  only  on  s  .  Sum- 
marizing, we  have  shown  that  cavitation  occurs  when  *]  is  greater  than  a  certain 
minimum  value,  v?    ,  that  this  minimum  value  depends  only  on  s  ,  and  that  cavi- 
tation is  produced  in  only  the  RG  combination  even  when  s  <  1  provided,  however, 

that  s  >  .05. 
o 

In  the  following  two  subsections  we   sketch  our  program  for  obtaining  the  de- 
pendence of    Y{,  s   /s1,    (ux  3)escA>0  and  (H     3/Hx)  on  (u     Q/bo)  for  various  values 
of  the  parameter  s    .     In  subsection    5J<2    it  is  assumed  that  s     >  1  or  that 
s    <  1  but  that    vj  -  1  >     (2y)/(y~1)   (s"J"-l)   (the  RG  combination);  in  subsection 
5^.3  it  is  assumed  that  sq  <  1  and  that    v^  -1  <     (2y)/(y-U  \(%  -1)   (the  RS  combin- 
ation).    When  s    <  1,  the  dependence  of  v^  ,  s./s^   (ux  e)esc/b0  and  (H     ^/H^)  on 

(u     „/b  )  for  all  u     „  >  0  is  obtained  by  combining  the  results  of  these  sub- 
y,0'   o  y,0 

sections.     The  results  are  summarized  in  the  graphs  of  Figures  15-l8  where  we 
have  plotted    v(  ,   s^,   (u     3)e8C/b0  and  (Hy^/H^  versus  uy^/bo  for  sQ  -   .0625, 
.1250,   .25,   .5,  1,  5,  10. 
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5.U2  The  program  for  the  RG  combination  s  >  1  or  s  <  1  but 


1  •*  £  &  Co"1-1'} 


t  = 


O 


Case  I:  >?<      Y{  . 
<           vcav 

From  RG,       and  RG7       of  (5.19)  we  have,   on  setting  u     9  and  u     „  equal  to 
zero  and  writing  s_  (/  0)  for  s„ 

rY(i-^  )s  -,5" 
(5.23)  H^I\V>\T)    "    h-D    —— -£] 

'_         1    +     n     YJ  J 


2  X 


V_    .  1   rStr(l-%av)(l+^  n  )    •]  ? 


5.2U)  ^jl.    -  -   ^ r_ J        X(s3/str,str). 


Recall  that  sA     «   >7  s     and  that  q         is  a  function  of  8 .     .     These  two  equations 
tr         '    o  ^cav  tr  ^ 

therefore  lead  to  relations  of  the  form 


RiH»  ST  '     b~  »  soJ  ' 

—  10  -1 


1,    2, 


which  do  not  depend  on  p  and  b  explicitly.  Now  we  assume  that  s  is  one  of  the 
values  mentioned  above;  our  program  for  establishing  the  dependence  of  vtt   s,/s,, 

K^esc^o  ■"*  (Hy,3/Hx)  on  ^,0^0  is  the  ^1°™%'- 


(1-1):         Choose   several  values  of   *\  <     V)  in  such  a  manner  that 

(%   »s+    )  "   1*1  s   >(  *\  s   )"   j  is  a  point  of  s  =  q"     and  at  the  same 
a  point  of  the  plotted  s(q,q.    :©)-curves    [see  Figure  6]. 

(a)     It  may  happen   (when  s     very  large  or  near  unity)  that  it  is  im- 
possible to  do  this  either  because  s.      is  too  large  to  be   shown  on 
the  graph  or  because  s.     <  1.0256  »   (.975)"  ,  that  is,   smaller  than 
the  smallest  value  of  8.      for  which  s(q,q.r,5)  is  plotted  on  the  graph. 
Under  these  circumstances  we  employ  the  linear  approximation  to 
s(q>cL   »5)  noted  in  Section  k   [cf.   (U.23)3  in  the  program  which  follows, 
The  dependence  of  p.         on  q.      may  be  obtained  from  Figure  8. 
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(1-2):  Compute  for  each  value  of  vi    chosen  in  1-1  the  value  of  the  expression 

for    ^(8o/8tr»8tr)  **  the  r ight-hand  side  of  (5.23).     Then  determine 

the  ratio  s  /s   from  the  plot  of  ^(s/s  js  )  versus  ^^(Figure  10). 

To  obtain  the  desired  accuracy  in  the  neighborhood  of  s,/s.  ^  0,  it 

5     tr 

proves  necessary  to  employ  suitable  approximate  expressions  for 
\\(s_/s   :s.  )  (cf.  (U.38)3  or  to  employ  a  graph  with  a  scales  that 
are  larger  than  those  employed  here.  The  same  is  true  of  X(s,/s.  ;s.  ) 
in  the  neighborhood  of  s_/s.   »  0  (see  Figure  11  and  (U.33)j] . 

(1-3):    Determine  from  Figure  11  the  value  of  \(s_/s  ;s     )  corresponding  to 

s_/s.   obtained  in  1-2. 
3  tr 

(I-U):    Determine  y\    and  s  /s   as  functions  of  u  _/b  .  To  obtain  the  de- 

pender.oe  of  vj  on  u  0/h  it  is  merely  necessary  to  compute  the  right- 
hand  side  of  (5.21)  as  a  function  of  >^  .  This  dependence  is  illustrated 
in  Figure  l5c. 

To  determine  the  dependence  of  s./s.   on  u  _/b  note  first  that  the 

5     tr  y»"     o 

values  of    v\    which  correspond  to  s_/s.     may  be  obtained  from  1-2.     The 

dependence  of  s  /s.      on  u     n/b     then  follows  from  the  relation  of  vi    to 
3     tr  y>*^     o  * 

u     _/b     Just  established. 

(l-5)t         Determine   (H     .,/H   )  as  a  function  of  u     Q/b   .     Note  first  that 

(H     ./H  )  ■    [(q~  -l)(l-q,s,)]   '     at  the  tail  of  the  rarefaction  wave 

and  hence  throughout  Region  (3).     In  Figure  7  we  have  plotted 

(H  /H  )  '    [(q^-Dd-q-.s.)]1/2  for  various  values  of  n       (=  s"  ). 
y     x  j  j  t>~  tr 

Since  s  /s       is  known,    (H     ,/H   )  may  be   obtained  from  this  graph.     When 
3     tr  y>*     *• 

s       is  very  large  or  close  to  unity,  Figure  7  no  longer  supplies  the 
necessary  information.     In  this  case  we  employ  the  linear  approximation 
to  s(q,q^    j 5)    [see   (U.23 )]]   in  conjunction  with  the  formula  for  (H     ~/H_J 
just  mentioned  and  compute   (H     ^/H  ).     The  results  are   summarized  in  the 
graphs  of  Figure  17b. 
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Case  II:  v?    >     v^. 

In  this  case  we  set  s0  equal  to  zero  in  (5.19)  RG-,       and  RG-,       and  employ 
the  second  relation  of  (5.20 ')  to  get  for   *]   >  Yi        : 


(5.25) 


i     J  r  f  U-H2)s„  ll/2         (        2    U/2  -r 


(5.26,      -g»°  .  i  J  JE-yas x(0i6tr; ♦  %1  (u  ,)eso 

o  I  U-     +  *\ 


u     0       u     ,.        , 

O  O 


'  p*  d-q     )(i+^ )    |i/2 

trv     Mcav/V     ^  '        I  . /n  , 

— -2— f        *(0'strK 


Here  the  ratio   (H     -,/H   )  at  the  tail  of  the  rarefaction  wave  is   given  by 
y,3     x 

the  expression     see   (5.19),  RG2;  take  s„  -  0 

(5.?7)  H^3/Hx     -    /q^d^)    ,  q3     -      q^. 

Assuming  s     is  one  of  the  values  employed  in  Case  I  above  we  proceed  in 
accordance  with  the  following  recipe: 

(II-1):       Obtain  the  value  of  X(C;s     )   for  each  s       '  \T     '     v\s0  from  Fi^1"6 


12  for  several  values  of      ^  >    v\ 


cav 


(II-2):       Obtain  for  each  q.      ■  s"     the  corresponding  q         «*  q~  from  Figure   8  and 
compute  the  ratio   (H     -,/H  )  by  means  of  (5.27). 

(II-3):  Obtain  the  value  of  ^(Ojs^)  from  Figure  12  for  each  str  -  ^  so  and 
then  compute  the  right-hand  side  of  (5.25)  for  each  vj  employed  above 
thereby  determining  the  dependence  of   (u     ,)-_./b     on  V]  • 

(Il-U)t       Utilize  the  information  in  II-l  -  II-3  and  compute  the  right-hand  side 

of  (5.26)  thereby  determining  the  dependence  of    >?    on  (u     n/bn)«     Graphs 

y>u  v 

of  the  relationship  between  v^  and  (u  0A>0)  for  several  values  of  the 
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parameter  s  appear  in  Figure  15c. 


(Il-5)s   Determine  the  dependence  of  (H  ,/H  )  on  (u  _/b  )  by  combining  the 

j>j    x     y>^>  o 

results  of  II-U  and  II-2.  Graphs  of  the  relationship  between  (H  ^/H  ) 
and  (u  _/b  )  for  several  values  of  the  parameter  e     appear  in  Figure 
17b. 

(II-6):   Determine  the  dependence  of  (u  ,)   /b  and  (u  ./b  )  on  (u  n/b  ) . 

x,j  esc  o      y»i     o      y*'-'  *» 

Fmploy  the  information  obtained  in  II-3  and  II-U.  Graphs  of 

(u  ,)   /b  and  (u  -/b  )  versus  u  _/b  are  raven  in  Figure  18. 
x,3'esc'  o      y,3  o'        y,0'  o 


5.U3  Tho  program  for  the  RS  combination  y\<   1  +  -^  (s~  -l)j  .0625  <  sq  <  1  . 
In  the  following  it  will  be  assumed  that  1  >  s  >  .C625.  Under  these  circumstances. 

(5.28)       y\<    i  +  |X-  (s"1-!)  <  v>    , 


cav 


(see  subsection  5.U1). 

Now  from  RS-       and  RS-       of  (U.2U)  we  have,  on  setting  u     „  and  u     „  equal 
->*x  _),y  x,£  y,^ 

to  zero, 

(5.29)  K(»3.>trntr)   -  r(str/si)1/2  <^> 

1  1         ! 

u     -       /p  \7  i  H    ,  v/sT  /  P«\Z  9" 

«*'  if  •(£)  I  ¥  I  *  t  (if )  '^/^w^- 

Assuming  s  fixed  ,  we  may  establish  the  dependence  of  yi  ,  (s,/s  )  and 

H  .,/H  on  (u  n/b   )  by  proceeding  in  accordance  with  the  following  prescrip- 
y>y    x     y,C  o'  J 


tion. 

(1-1):  Choose    ^>  j.  <yj<  1  +  ^-  (s"1-!),   such  that  the  point   (q^s^  =   (q^  v?sq) 


2y  /  -i 

— T    (S   - 

Y-l   o 
is  a  point  of  intersection  of  an  s(q,q.  :5)-curve  and  the  curve  whose 


*  In  our  computations  we  chose  s   to  be  .5,  .250,  .125  and  .0625. 
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parametric  equations  are  given  in  (5.13)  and  (5.1U).  Each  such  point 
may  be  obtained  graphically  from  Figure  13  or  an  enlarged  version  of  it, 

(1-2):    Obtain  o„   and  a    corresponding  to  (q,,s, )  from  Figure  13. 
^cr      csv  -L  j» 

(1-3):    For  each  s  determine  (P1-P£J)/P0  from  Figure  U. 

(I-U):  Compute  the  right-hand  side  of   (5.29)  and  then  find  s_/s     from  the 

plot  of    \\(s  /£>js,    )  versus  s^/s,.     (These  curves  are  not  given  here; 
they  may  be  obtained  from  the     ^(s-./s,;s.    )-curves  of  Figure  10  by 
making  use  of  (U.20).) 

(1-5):         Determine  the  value  of  X.(s^,s,;s     )  from  the  \(s_,s1;str)  versus  s.,/s.j- 
curves.     (These  curves  are  not  given  here;  they  may  be  obtained  from 
the  X(s-/s     ;s     ) -curves  of  Figure  11  by  making  use  of  (U.20).) 

(1-6):         For  each  sn    (»  >7  s    )  determine   H     ,/H     from  Figure  5« 

J.  >    o  y  ,a    x 

(1-7):         Compute  the  right-hand  side  of  (5.30),  for  each  value  of    yj  employed 
in  the  previous  steps  and  obtain  in  this  manner  the  dependence  of  Yj 

on  u    Q/b0,  1  <  v]<    Y}crit  -  1  +  rjr-  (Sq1-1^     In  order  t0  keep  the 
diagram  uncluttered,  the  >j -versus   (u     _/b) -curves  in  this  region 

have  been  omitted  in  Figure  l5c.     However,   these  curves  may  be  obtained 

from  Figure  15b,  where  we  have  plotted  the  corresponding    vi    versus 

u     »/a     =  u    n/b    Vs~"  -curves. 
y>0     o        y,0'    o  r    o 

(1-8):    Compute  H  ^/K%   =  U^-lKl-q^M.  Note  that  the  value  of  sy^  is 
known  from  I-U.  Compute  s.,  ■  (s_/s,)s,  and  determine  q^  from  the 
appropriate  s(q,q.  ;5)-curve  in  Figure  13  [see  1-2]. 

(1-9):    Determine  the  dependence  of  (s./s-^  on  (u  Q/b0)  by  plotting  the  values 
of  (s  /s  )  obtained  in  I-U  versus  the  values  (u  o^o^  obtained  in  I-7 
(see  Figure  16] • 
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(1-10):       Determine  the  dependence  of  (H     -i/fL)  on  u    0/b  j 

0  i  'uy  o^o^  -  ^uy  O^crit^bo'  ^  combinin?  the  results  of  1-8  and  1-9 
[see  Fi  gure  17bJ  . 

5.5     Approximate  expressions  for  vj  ,    sjs,    and  (H     .,/H  ) 

5.5l     Approximate  expressions  for  v\  ,    (H     -j/H-J  when  (u     n/b   )  is  very 
y ,  i     x     y,u     o     

large.     Let  us  choose  u     _/b.    so  large  that  for  the  corresponding  -q   we  have 

(5.31;  ^  »     ^  >     1j       whenever  s     >   .05 

v  '  l  cav  '  o  — 

or 

(5.32)  yj   »  Y,orU  -  1  ♦  ^  (8^-1;  >     r\  cay  >  1,     whenever  so  <   .05. 

In  both  cases,  as  we  have   shown  in  subsection  5.U1,  the  resolution  of  the   shear 
flow  discontinuity  is  effected  by  an  RG  combination.     For  fixed  s     we  may  assume* 
since     vi    is  very  large  compared  with  unity,  that 

(5.33)  str    -     y\   so  »  1 
and  consequently  that    (cf .   (U.28JT] 

tr 

Since     *j    exceeds     V)        ,  the  quantity  s.  vanishes  at  the  tail  of  the  rarefaction 
wave  and  hence,   as  already  noted  in  (5.27), 

1/2 


(Hy,3/Hx)     -      l^cavKav] 


Employing  (5.3U)  we  then  conclude  that 

0+l  1/2  1/2      0+1  1/2        1/2 

(5.35)        (Hyj3/Hx)  a  (2ji)   sj/2  -  <2£)   (nso) 

Neglecting  lower  order  terms  we  have  from  (5.20),  (5.21)  and  (5.35) 
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b  =       H       v   x,3  esc 

o  x 


(u     ,)  ,      2       1/2  1/9 


Consequently,  we  may  conclude,   since  &  *  y/(2—f)t  u-  «*   (y-1)/(y+1),   that 
(5.36)  1  «0   «4   J  U*Y)1/2   (   ^  ) 


(5.37)  ^3  ,    (UY)1A  (  ^0  ,1/2 

X  o 


k  -*\                    (ux,3}esc  v  i         /  uy,o  ,1/2 

(5.38)  _ g  lA    (  -£-  J 

o  (1+y)  ° 


Moreover,  since 

H 


u 


^  (u  ,) 


y>3     H    x,3  esc 


at  the  tail  of  the  rarefaction  wave  we  have, 


(5.38-)  ill    „    !zi0 

o  o 


When  the  adiabatic  exponent  y  is  taken  to  be  5/3  we  find 

(5.36-)  ^6     $   1.36  (   -|l°   ), 

o 


(5.37M  iil  s  1.27  (  ^  )V2  , 

X  o 

(5.38,)  (V)eSC   a   .787  (  V  )1/2  . 


bo 


The  relation  (5.37')  replaces  the  incorrect  result  for  a  monctonic  gas 
given  in   [7],    (see  pp,  U9-5CQ . 
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5.52     Approximate  expressions  for  vo  ,     s-j/s+     gJld   (H     ^/H^)  when 


(u     „/b   )  is  small. 
y,0'    o'  


Case  I:      s   >  1  (the  RG  combination) 


When  s  >  1,  the  resolution  of  the  initial  discontinuity  is  effected  by 
the  RS  combination.  Assuming  *]  -1  is  small  we  have  from  (5.23)  and  (5.2U), 
on  neglecting  terms  involving  higher  powers  of  (  Y|  -l), 


(5.39) 


(5.1*0) 


Since 


^VW     '     (V-L>(s<//2 


rl/2 


JLl£     S     i       B   (i.q       )  X(s./8.    js+    ), 

b  y   I      °         cav'J  3'    tr*  tr 


(5.U1) 


,8tr     "     X 


when    V)  -1  and   (u     n/b   )  are  small    [see  Figures  15  and  l6] ,  the  relations   (5.39) 
and  (5. U0)  become,   in  virtue  of  (U.UO)  and   (U.I4I), 


(5.U2) 


(5.U3) 


1-  (s3/Btr)]    a    (n-D(s0)1/2 

U  p2s     (1-CL.  J-,1/2,-  -r 


1/2 


Now 


l-qtr  3  l-(^  s0)-L|_(vl-l)s0  ♦   (s0-l)]  / 


Sd' 


Consequently 


(5.U3-) 


u 


(  Ml  )  *  |" 


2{(vl-l)80+(8-l)}-,l/2 


T  &«>*>¥* 


Equations  (5.U2)  and  (5.1*3')  when  solved  for  (^   _i)  and  l-(s  /str)  yield  the 
following  results 
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(5.UU) 


Y(-l 


,      s   -1 


2Ys 


1/2 


-1+^1  + 


(ae-D4 


o 


(5.U5) 


hw]-  ^ 


s  -1 


-1  +J1 


(8    -1)^  O 


When  s     >  1  and  (u     _/b  )  is  sufficiently  small  we  get 
o  y,0'   o  -v  b 


(5.U6) 


^v1' 


(5.U7) 


VStr  "  1  "  51^ 


=ir(x-)  ■ 


s      >  1. 
o 


On  the  other  hand,  when  s  is  unity  we  have 


(5.U8) 


Y|   S  1  +   a}l/2  (  Vo   , 


(5.U9) 


u 


we 


Now  since  s      ~    s       we  must  have  q.  ~  q^   .     Employing  this  fact  and   (U.U1), 
obtain  for   (H     _/H  )    (cf.  RG„    (5.19)1 

r,  nV2r  l1'2 


Comparing  this  expression  with  that  for  (u     fl/b   )  in  (5.U3)  and  noting  that 
So(l-qtr)  "    |Bo^8tr"1,/Vl   So]  ~  ^tr"1^  we  have  finally 


(5.30) 


(Hy,3/fV    *    (uy,0/V> 


(cf.    t7],  p.   50). 


-  72  - 


Case  II:     s   <  1  (the  RS  combination) 
o 

When  8  <  1,  the  resolution  is  effected  by  an  RS  combination.  Now  when 
o 


"ri   -i  is  sufficiently  small  £  and  S„  of  (5.2)  become 


(5.51) 


pl"po 


m=L 


^-1*0 


/H   .\2     2(p  -p  )  2U"0 

(5.52)  \^f±)       *  _±-°-  (l-so)«— -^  (T-i),  yj-i^o 


Moreover,  when 


(5.53) 

and 
(5.5U) 


>ls       <     1 
1    o 


8       ■       S, 


the  expressions  for  K(s  /S-js.  )  and  A.(s.,/s  js.  )  may  be  approximated  as  follows: 

j     l  xr        j     j.  xr 

u 


(5.55) 


l\(s3,si;str)  -  j  (q/V1'^  ?d^^^V\r,1/2E  "  (VSi)] 


^(83,8^8^) 


(l-q)/(l-qcav) 
1  -  qs 


(f)"5d(f-) 

Sl     Sl 


(5.56) 


P^JH*>] 


Therefore  when  vj  -1  is  sufficiently  small  we  have  from  (5.29),  (5.30)  and 
(5.51),  (5.52) 


(5.57) 


|i  -  (ys±)]   -   (siqi)-l/2(yi  -  1) 
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(5.59) 


V,0 


'2(l-soH^-l) 


"  o 

Y 


i-q, 


l-q^i 


_l-(s3/si)] 


where  [cf.  (5.13)J 


(5.59) 


(^MH^)2-  [(^1)2^(^i)(Hy,1/Hx/  +  0((H^i/Hx)U}] 


1/2 


Neglecting  the  term  0  |(H  , /H  )  f  ^d  assuming,  in  addition,  that  the  ratio  of 
the  second  and  the  first  term  in  the  braces  is  small  compared  with  unity,  we  con- 
clude that 


(5.60) 


or 


(5.60') 


h    a  (i-i^-i)],   r,- 


1  -\  -  l^n-1^     n-1  "°- 


1  -  0 


Substituting  this  expression  for  (1-q, )  and  the  expression  for  l-(s-/s.)  given 
in  (5.57)  into  (5.5b)  we  get 


(5.61) 


v.  [^](V1)?-H^r^'7]: 


1/2 
As  a  result  we  have,  on  solving  for  (>o  -1)  '    , 

4 

i-  Y(1-so) 


(5.62)  ty-iy 


i  -  a  -  U(I) 


Y>2"         1 
Y(l-so) 


1 


(4^) 

o 


Assuming,   the   second  term  in  the  brace  is  small  compared  with  unity,  we  get 


<5-«>        n^-zuV  ft2)2- 


U 


y,o 


<-*>  0 


and,   in  addition, 


-  7U  - 


J  2s  '    (1-s    )      \      o    /  o 


rare- 


Now  the  ratio   (H     3/H   )   -  Kq^-lJd-q^J ip2  at  the   tail  of  the 

faction  wave  is  effectively   (H  i/Hx)>   since  q_  ^  q.^  and  s    ^  s.^     Consequently, 
we  have   from  (5.63)   and   (5.52) 

(see    [7],  P.   50). 
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6.       Interpretation  of  the  results 
6.1     Introductory  remarks 

On  the  basis  of  the  results  of  the  preceding  section  we  are   in  a  posi- 
tion to  obtain  a  complete  description  of  the  motion.     More  explicitly,  by  choosing 
successively  greater  values  of  time,  t,  we  could  obtain  a  picture  of  how  the 
various  disturbances  propagate  away  from  the  initial  discontinuity    [see  for  ex- 
ample Figure  T\,   and  by  varying  the  parameters  which  characterize  the  state   in 
front  of  the  shocks  we  could  establish  quantitatively  the  dependence  of  the  solu- 
tion in  each  of  the   flow  regions  on  the  state  of  flow  in  the  undisturbed  medium. 
To  carry  out  such  a  program  adequately  would,  however,  be  beyond  the  scope  of 
this  report.     We  shall  instead  limit  ourselves  to  discussing  certain  general 
features  of  the  motion.     Attention  will  be  focused  exclusively  on  the  character 
of  the  motion  in  the  region  behind  the   shock  and  in  the  region  at  the  tail  of 
the  simple  wave. 

Let  us  recall  some  of  the  notations  and  conventions  employed  in  the 
preceding  sections.     First,  in  view  of  the  symmetry  of  the  motion,  there  is  no 
loss  in  generality  in  limiting  our  considerations  to  the  motion  along  the  positive 
x-axis.     In  Figure  3a  we  have  depicted  the  basic  flow  regions  in  the  positive  half 
of  the  (x,t)-plane.     These  basic  flow  regions  are:  Region  (0),  the  region  of  con- 
stant state  in  front  of  the  (forward-facing)   shock;  Region  (1),   the  region  of  con- 
stant state  behind  the  shock;  Region  (2),  the  flow  region  of  the  slow,  centered, 
forward-facing  rarefaction  wave;   Region  (3),  the  region  of  constant  state  following 
the  tail  of  the  forward-facing  simple  wave.     The  state  in  Region  j,   j  »  0,  1,  2, 
3  is  characterized  by  the  quantities  p.,  b       -  ('jHxP1   )  >   sj     "^ai/'bj^  *  ^Pj/^  * 

(Hy,j/Hx)j  "x,^  (or  **, j/*a>  «*  y^j  (or  vyv-  Here'  in  Region  (3)'  PJ 

is  the  density;  b  the  Alfven  disturbance  speed;  s  the  square  ratio  of  the  sound 
speed  a  to  the  Alfven  disturbance  speed  b  ;  H   .  and  u   .  the  components  of  magnetic 

j  j     y » j    y  >  3 

field  and  fluid  velocity  respectively,  transverse  to  the  x-axis;  and  ux  .  the 
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longitudinal  component  of  the  fluid  velocity.  In  addition,  H  is  the  fixed 

longitudinal  magnetic  field,   [i  the  magnetic  specific  inductive  capacity,  and  y 

■a- 
the  adiabatic  exponent   ,     In  terms  of  these  notations  our  chief  objective  may 

be  expressed  as  follows:     To  find  in  Regions  (1)  and  (3)  qualitative  information 
about  the  dependence   of  the   solution  on  (i)   the  ratio  u     0/aQ  for  fixed  so,  and 

on  (ii)   8  («  YP0/mh£)  for  fixed  uy  0/ao  °r  f°r  fl*ed  ^O^o*   accordine  as  80  is 
made  to  vary  by  varying  H    or  pQ.     It  will  be  recalled  that  u     Q  is  half  the 
magnitude  of  the  initial  discontinuity  in  the  transverse  velocity.     It  should 
be  emphasized  that  the  corresponding  information  could  easily  be  obtained  for 
the   simple  wave  region,  Region  (2),  but  for  the  sake  of  brevity  we  shall  not 
do  so. 


6.2     Some   general  features  of  the  solution  in  Region  (1),  behind  the  shock 
The  dependence  of  the   state  in  Region  (l),  behind  the  shockjOn  s     and 
u    n/a     or  any  equivalent  set  of  parameters  may  be  obtained,   as  we  have  already 
observed  in  Section  5,  from  the  knowledge  of  how  the  shock  strength    vj  -  Pt/P0 
that  is  the  ratio  of  the  pressure  p,  behind  to  p     in  front  of  the  shock,  depends 
on  these  parameters.     As  is  perhaps  to  be  expected,  the  dependence  of   vj   on  these 
parameters  may  differ  according  as  the  resolution  is  effected  by  a  gas  shock  or 
a  switch-on  shock  or,   stated  analytically,  according  as 

(6.1)  so    >     1         or         so  <  1        but      V.     v?crit 
or 

(6.2)  so    <     1,       but    ■      V|<    Y^   . 

It  will  be  recalled  that 


*  It  should  be  recalled  that  all  graphical  results  are  based  on  the  choice  of 
Y  ■  5/3  for  the  adiabatic  exponent. 
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<6'3)  7crit     '     1  +  |?T     (^1-1^  so     <     l- 


is  that  value  of  the    shock  strength  above  which  the  svitch-on  shock  reduces  to 
a  pure  gas  shock.     In  the  following  we  shall  treat  these  two  cases  separately. 

6.21     Some  general  features  of  the  solution  in  Region   (1)   behind  the 
pure  gas  shock.     Here,   in  accordance  with  (6.1),  we   assume  that  s     >  1  or  s     <  1 
but     "V{  >  Vj  Consider  Figure  l5a>where  we  have  plotted  p,/p     ■   v?  ■ 

=  Y|  [u    n/a  ,s  "J  versus  u     -/a     for  several  values  of  the  parameter  s   .     For 
values  of  s     <  1  we  have  plotted  only  those  portions  of  the     v?    versus  u     n/a  - 
curves  which  correspond  to  values  of  v^  which  are   greater  than       vi  From 

the  figure  we  conclude  that  for  fixed  s     the   shock  strength    v?    increases  mono- 
tonically  with  u     _/a  •     In  addition,  as  we  have   shown  in  Section  5,    (5.36),   (5.U6) 
and  (5.63), 

1/2  2 

<6-*>        -i  ■  *  ♦  j£nr  (¥■)  •        ^  -  °.       so  >  i 

o  \      O     '  0 


(6.U0  *?  -     1    +(J)1/2^V  ^    ■*    0, 

V         O      /  O 


(6.5)  v(  -    J    (^C)1/2    i£    ,  ^    -  oo,  sQ  >  0. 


Since   (p,-p   )/p  ,  u     ,/a  .and  U  /a     increase  with  increasing  vj      [cf.  G, ,  G-,   G,  , 
v   1     o       o'     x,l'    o>  x     o  I      *-         1       3       U 

(5.5)3  we  mav  conclude,  moreover,   that  these  quantities  increase  with  increasing 

u     0/a  ,  assuming  s     is  fixed.     Both  u     ,/a     and  U  /&     approach  infinity  with 
y,0     o  o  x»l     oxo 

2     2 
u     ./a    while   (p-.'P   )/p     approaches  the  constant  value   (l-u.  )u.     »  3   (when  y  ■  5/3) 
y  ,u     o  ±     o       o 

as  u     n/a     approaches  infinity.     The  transverse  velocity  is  continuous  across  a 

gas  shock    [see  G^     (5.5)]  while  the  transverse  magnetic  field  vanishes  on  both 
sides  of  the  shock  front. 
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Thus,  if  the  values  of  the   sound  and  the   Alfven  disturuance  speeds 
in  front  of  the  pure  gas  shock  are  held  fixed,  the   values  of  the   shock  strength, 
the  excess  density  ratio,   the  speed  of  the   shock  front^and  of  the  components  of 
fluid  velocity  behind  the  shock  vary  in  the   same   sense  as  the  magnitude  2u 
of  the   initial   shear  flow  discontinuity . 

Now  let  us  fix  u     _,   the  pressure  p    ,and  the  sound  speed  in  front  of 

y,0  o »  r 

the  gas  shock  and  vary  H  .  Assume,  in  particular,  that 

(6.6)  Hx  >  Hx  . 

Then  the  ratio  u     ./a     is  fixed  but  the  quantity  s     varies,  i.e., 
y,0     o  *     o  *  * 

(6.7)  s^     -     YP0/^    <     sQ     -     YP0/HHX   • 

If  now  f  ■  f  (u     -/a   ),s       is  any  function  of  (u     _/a   )   and  s   ,   let  f  denote   the 
|_  y»0     o   '   oj  J  v  y,0'    o'  o* 

value  f     (u     n/a   ),s     .     Then  from  Figure  15a  and  from  G,,  G.,  and  G^,    (5.5)  we 


conclude  that 


^o  Ko 

(6.3) 

u     ,   >  u     ., ,       U     .,     >    U     ,    • 
x,l         x,l'         x,l  x,l 

Thus,  if  the  magnitude  of  the  initial  shear  flow  discontinuity  and  the  values 
of  the  pressure  and  the  density  in  front  of  the  gas  shock  are  held  fixed, then  the 
values  of  the  shock  strength,  the  excess  density  ratio,  the  shock  speed  and 
the  longitudinal  component  of  particle  velocity  behind  the  shock  vary  in  the  same 
sense  as  the (fixed)  longitudinal  component  H  of  the  magnetic  field,  and  hence  in 
the  opposite  sense  to  the  square  ratio  s  of  the  sound  speed  to  the  Alfven  dis- 
turbance speed  in  front  of  the  shock. 

In  Figure  15c  we  have  plotted  >(  ■  *\  u  0A>o>so  versus  u  Q/bo  for 
several  values  of  s  .  As  in  Figure  15a,  the  portions  of  the  curves  corresponding 


It  is  assumed,  of  course,  that  the  Y[    -  ■*}(£  )  and  ^  -  >v?  (s  )  are  both 

greater  than  critical  values  corresponding  to  &  and  s  . 

o     o 
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to  s  <  1  and  lying  below  the  critical  value  have  been  omitted.  Employing  this 
figure  and  G. ,  (5.5)  and  proceeding  in  the  above  manner  we  conclude  that  if  the 
magnitude  of  the  initial  shear  flow  discontinuity  and  the  values  of  the  longitu- 
dinal magnetic  field  and  the  density  in  front  of  the  gas  shock  are  held  fixed,  then 
the  shock  strength  and  the  excess  density  ratio  increase  as  the  pressu re  in  front 
of  the  shock  decreases  and  hence  as  the  square  ratio  s  of  the  sound  speed  to 


the  Alfven  disturbance  speed  in  front  of  the  shock  decreases. 

Unlike  the  previous  case,  the  qualitative  dependence  of  u  ,  and  U  , 

X  y  JL  X  «  J- 

on  p  and  u  _/b  is  not  immediately  revealed  by  G,   and  G^,  (5.5)»  This  de- 
pendence could  of  couise  be  established  graphically,  but  we  shall  not  do  so  here. 
It  was  stated  in  Section  1  that  the  motion  of  the  medium  differs  radi- 
cally from  the  corresponding  steady- state  gas  dynamical  shear  flow  motion  even 
when  H  is  small, if  the  magnitude  of  the  initial  shear  flow  discontinuity  is  large 

enough.  This  statement  follows  from  the  fact  that  the  value,  2(u  n)   ,  of  the 

y,u  cav 

initial  discontinuity  at  which  cavitation  occurs  remains  finite  as  H  approaches 

zero  and  hence  as  s  vanishes;  all  other  parameters  are  assumed  fixed.  Actually, 

(u  „)   /a  approaches  the  value  3»67  as  s  approaches  infinity,  a  fact  which 
v  y,0'cav'  o  PK  o  ^r  *" 

may  be  inferred  from  an  inspection  of  Figure  l5&.  Clearly  then,  in  contrast  to 

the  corresponding  gas  dynamical  situation  where  the  motion  is  steady  and  purely 

transverse,  the  resolution  of  the   shear  flow  discontinuity  involves  a  longitudinal 

motion  which  results  in  cavitation  even  for  very  small  H  if  u  ~/a^  is  somewhat 

J  x    y,0'  o 

larger  than  the  value  3»67« 

The  radical  difference  between  the  gas  dynamical  motion,  which  may  be 
regarded  as  a  hydromagnetic  motion  in  the  absence  of  magnetic  fields,  and  the  hydro- 
magnetic  motion  considered  here  could  perhaps  have  been  anticipated  from  our  dis- 
cussion of  the  shear  flow  layer  in  subsection  2.(j.  We  showed  there  that  the  force 

2/OU   \/d  u   \     ?_ 

per  unit  volume  F,  within  the  layer,   is  given  by  the  expression  -y^Z(-^-)( h)  t  x  > 

ox 
where  t,  the  time,   is  sufficiently  small.     Clearly,   no  matter  how  small  H    is   tsken, 
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it  is  possible  to  make  T  as  large  as  one  pleases  by  making  the  3u/9x  large 
enough.     This  may  be  accomplished  by  making  the  difference,   2u     „.  of  u     at  the 

y»o*         y 

rigfict  and  the  left  side  of  the    shear  flow  layer  large  enough.     Thus,  by  making 
u     q  large  enough  one  might  expect  to  produce  a  large-scale  longitudinal  motion 
and  perhaps  even  cavitation, as  long  as  H     does  not  vanish  altogether. 


6.22  Some  general  features  of  the  motion  in  Region   (1)  behind  the  switch- 

on  shock.     In  accordance  with  (6.1)  we  assume  here  that  s     <  1  and    "0  <     >1       ,  .• 
o  <    _      i  cnt 

In  Figure  l5b  we  have  plotted    Y|  -  ^    (u       /a   ),s         versus   (u       /a   )   for  several 
values  of  parameter  s     less  than  unity  and  for  (u     _/a  )  in  the  range' 

0  -uy,0/ao  1  fry^eril/V     Here>   (uy,0}crit  is'  for  fijced  V  the  value  of 
u     Q  when     X\  ■   v?  t     When  u         is  small  we  have    (cf.   (5.63)] 

Vs  n        /u     -x2 


<6'9>  vj«     i     +     ?(T^1 


o     f  y>° 


o-     \   ao 


From  Figure  15b  and  S,,  S,   and  S.  ,  (5.2)  we  deduce  the  following  facts. 
If  the  values  of  the  sound  and  the  Alfven  disturbance  speeds  in  front  of  the  switch- 
on  shock  are  held  fixed,  then  the  shock  strength,  the  excess  density  ratio,  and  the 
longitudinal  component  of  the  particle  velocity  behind  the  shock  are  monotonically 
increasing  functions  of  the  magnitude  of  the  initial  discontinuity  in  the  transverse 
velocity,  wnereas  the  speed  of  the  shock  front  decreases  monotonically  as  this 
magnitude  increases.  Furthermore,  if  the  magnitude  of  the  initial  discontinuity, 
2u.  0,  and  the  values  of  the  pressure  and  the  density  in  front  of  the  switch-on 
shock  are  held  fixed,  then  the  shock  strength,  the  excess  density  ratio  ,  the 
longitudinal  component  of  fluid  velocity  behind  the  shock,  and  the  speed  of  the 
shock  front  increase  with  decreasing  H  and  hence  with  increasing  s  . 

It  should  be  observed  that  in  a  pure  gas  shock  the  variation  of  >o  , 

(p.-p  )/p  ,  u  ,  ,and  U  with  H  for  fixed  u  _  and  a„  is  opposite  to  the  variation 
Kl  o"   o*  x5l'     x      x  y,0     o 

*  The  dependence  of  (pn-p  )/p  on  s  was  established  graphically. 
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of  the  corresponding  quantities  in  a  switch-on  shock    [see  subsection  6.21~l.     In 

addition,  whereas  u       -  u         and  H     ..    •  0  behind  the  pure   gas  shock,  in  the  pre- 

y  >-*■     y>u  y»-«- 

sent  case  u  -  and  H  ,  differ  from  the  corresponding  values  in  front  of  the 
shock  and  moreover  depend  on  H  .  The  nature  of  this  dependence  could  be  ob- 
tained from  S2  and  S^     ,  (5.2),  and  Figures  U,  5,  and  l5bj  however  we  shall  not 

include  these  considerations  here. 

If  we  were  to  superimpose  Figures  I5aand  15b,  the  family  of  curves 

corresponding  to  values  of  s     >  1  would  intersect  those  corresponding  to  values 

of  s     <  1.     It  is  therefore  not  possible  to  make  any  general  statement  regarding 

the  behavior  of  the  entire  family  of  curves  which  would  be  valid  in  the  common 

range  of  u     n/a  ,      It  might  appear  strange  that  the  pressure  p-,/p     is  smaller 

behind  switch-on  shocks  than  behind  gas  shocks  when  u     ./a     is  sufficiently  small. 

y,0'    o 

This  fact  could  perhaps  be  explained  by  recalling  that  the  total  pressure  behind 
a  switch-on  shock  includes,   in  addition  to  the  usual  pressure  term,   a  term  in- 
volving the  magnetic  pressure,   \ih     ../2,and  that  the  contribution  of  this  term 
to  the  total  pressure  behind  the   shock  is  not  taken  into  account  in  our  definition 
of  shock  strength. 


6.3     Some   general  features  of  the  motion  in  Region  (3)  between  the 

rarefaction  waves       Turning  now  to  the  discussion  of  the   state  between  the  tails 

of  the    simple  waves,  let  us  first  consider  Figure  19.     Here,  we  have  plotted 

p„/p     ■  b^/s     versus  u     ^/b     for  several  values  of  the  parameter  s   •     The  curves 
3     o         y   o  y,0'    o  ^  o 

of  this  figure  may  be  ohtained  from  those  of  Figure  16  by  using  the  relation 
p_/p     =    ^(p^/P-,)   =    ^(s./s   ).     In  contrast  to  the  curves  of  Figure  16,   however, 
the  curves  of  Figure  19  corresponding  to  different  values  of  s     do  not  intersect. 
It  is  easily  shown  that  the  p./p  -versus  u     0/b  -curves  enjoy  the  same  property 


and 


moreover  vary  with  u  _/b  in  the  same  way  as  do  the  p^/p  -curves.  We  con- 


4t 

An  approximate  expression  for  p./p  when  u  _  is  small  may  be  obtained  from 
the  results  of  subsections  5.51  and  5.52. 
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elude  therefore  that  if  the  magnitude  of  the  shear  flow  discontinuity,  2u  _, 
and  the  values  of  the  longitudinal  magnetic  field  and  the  density  in  front  of 
the  shocks  are  held  fixed,  then  the  pressure  p  in  front  of  the  shocks  and  the 


ratios  p  /p  ,  P-j/P0  vary  in  the  same  sense.  It  will  be  recalled  that  p,  and  p. 
are  the  pressure  and  density  in  the  Region  (3)  behind  the  rarefaction  waves 
while  p  and  p  are  the  corresponding  quantities  in  the  undisturbed  region  in 


front  of  the  advancing  shocks* 

Similarly,  from  an  inspection  of  Figures  l?b  and  19  we  conclude  that 
if  the  magnitude  of  the  shear  flow  discontinuity  and  the  values  of  the  longitu- 
dinal magnetic  field  and  the  density  in  front  of  the  shocks  are  held  fixed,  then 
the  pressure  p  in  front  of  shocks  and  the  magnitude  of  the  longitudinal"  escape 
velocities"  and  of  the  transverse  velocities  at  the  tails  of  the  simple  waves 
vary  in  the  opposite  senses  whereas  the  pressure  in  front  of  the  shocks  and  the 
transverse  component  H  ,  of  the  magnetic  field  in  the  region  between  the  rare- 
faction  waves  vary  in  the  same  sense . 

From  Figure  1? a, where  we  have  plotted  H  -,/H  versus  u  „/a  it  is  an 
E      '  e  y,3  x       y,Cr  o» 

easy  matter  to  deduce  the  following  result.  If  the  magnitude  of  the  initial 
discontinuity  2u  _  and  the  values  of  the  pressure  and  the  density  in  front  of 
the  shock  are  held  fixed,  then  the  transverse  magnetic  field  H  ,  in  the  region 
between  the  rarefaction  waves  is  a  non-decreasin?  function  of  the  (fixed)  longi- 
tudinal magnetic  field  H  . 

When  u  n  is  small  we  know  from  (5.50)  and  (5.65)  that 

(6.10)  H  ,/K   *  u  _/b  u  .  ->  0 

y>3'  x     y,0'  o  y,0 

or}  equivalently}  that 

(6.10.)  Hyj3    *    (PoA01/2u7)0  uy>0  -    o. 

On  the  other  hand  when  u  n  is  large  we  know  from  (5.37)  that 

y,u 
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(6.31)  Hy;3/Hx     •     (Ur)a/U(uy;0/bo)1/2,  uy>0  ->  oo 

or,  equivalently,  that 

(6J1-)  Hy>3     S      [(l+Y)po/qiA  (H^)1'2,         uyj 


0  "*   °°' 


Thus,   H     .,  is  independent  of  H    when  u     _  is  small  and  increases  as  the  square 

y,3    tZ x    y,0    3 

root  of  H  and  the  fourth  root  density  when  u  _  is  large. 
x  * y,0  a- 

It  should  be  emphasized  that  care  must  be  exercised  in  applying  (6.11') 
to  a  given  physical  situation.  The  reason  for  this  is  the  following:  Equation 
(6.11"')  is  valid  only  for  sufficiently  large  values  of  u  n  and,  in  particular, 
only  for  those  values  of  u  n  which  greatly  exceed  the  minimum  value  of  this 
quantity  necessary  to  produce  cavitation.  Even  if  it  be  granted  that  it  makes 
sense  to  use  the  fluid  dynamical  notions  of  fluid  velocity,  density  and  the  like, 
the  objection  to  assuming  that  the  conductivity  of  the  fluid  is  infinite  must  be 
met. 

As  we  have  already  emphasized  in  the  introductory  section  of  this  work, 
(6.11')  implies  that  arbitrarily  large  transverse  magnetic  fields  may  be  created 
in  the  region  behind  the  rarefaction  waves,  provided  the  fixed  longitudinal  magne- 
tic field  component  does  not  vanish  altogether  and  provided  the  magnitude  of  the 
initial  shear  flow  discontinuity  is  made  sufficiently  large.  Associated  with 
this  is  the  fact  that  the  motion  which  resolves  the  initial  discontinuity  may, 
under  the  same  circumstances,  be  made  to  differ  radically  from  the  corresponding 
gas  dynamical  motion.  In  fact  shocks  of  arbitrarily  large  strength  may  be  pro- 
duced since  (cf.  (6.5)3  the  expression  ratio  p, /p  of  the  pressure  p,  behind  the 
(gas)  shocks  to  the  pressure  p  in  front  is  given  by  the  following  expression: 
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See  Fig.  3  on  following  page 

FIG.  3.  Waveforms  of  a  resolution  consisting  of  fast  gas 
shocks  followed  by  slow,  centered  rarefaction  waves.  The 
fast  gas  shocks  must  be  replaced  by  switch-on  shocks  when 
the  Alfven  disturbance  speed  is  supersonic  ahead  of  the 
shocks  and  the  magnitude  of  the  initial  discontinuity,  2u   , 
is  below  a  certain  critical  value  which  depends  upon  the 
ratio  between  the  sound  speed  and  Alfven  disturbance  speed. 
In  this  case  the  ratio  (H/H  )  [see  Figure  (c)_]  of  the  trans- 
verse to  the  longitudinal  component  of  the  magnetic  field 
is  positive  in  the  region  behind  the  shocks.  It  is  assumed 
in  the  above  diagrams  that  the  magnitude  of  the  initial  dis- 
continuity is  below  the  minimum  value  required  to  produce 

cavitation.  When  2u  n   exceeds  this  minimum  value,  p.,  the 

y,0  Y 

pressure  in  Region  (3)»  vanishes.  In  addition,  the  longi- 
tudinal component,  u  ,  and  the  transverse  component,  u , 

x  y 

of  the  particle  velocity  no  longer  vanish  at  the  "  tails" 
of  the  simple  waves  -  that  is,  at  the  inner  boundary  lines 
of  the  simple  waves  [see  Figures  (c)  and  (e)J}  both  velo- 
cities are  positive  in  the  region  x  >  0  and  negative  in 
the  region  x  <  0. 
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FIG.  6.     Graph  of  the  relationship  between  s  ■  s(q,q.    ;5)  and 
q  in  a  slow  simple  wave.  The  ordinate,   s,   is  the   square  ratio 
of  the   sound  speed  to  the  Alfven  disturbance   speed,  while  the 
abscissa,   q,   is  the   square  ratio  of  the  slow  simple  wave  dis- 
turbance speed  and  the   sound  speed.     The  function  s(q,q,    j5) 
is,   for  each  q^    ,  the   solution  of  the  basic  ordinary  linear 
differential  equation   (U.S).      JAs  usual  it  is  assumed  that 
Y,  the  adiabatic  exponent,  is  5/3»   hence  0  =   5.Q  After  inte- 
gration (U.3)  leads  to  the  expression  given  in  (U.13).  Note 
that  when  q  ■  q^   ,  then  s  ■  s       -  q"     .     Thus,   (*Lr>s+r)  is 
a  point  of  the  arc  s  ■  q"   .  The  various  s-versus-q-curves 
correspond  to  the  following  choices  for  the  value  of  the 
parameter  q^:  q^  =■   .975,    .95,   .925,    .9,    ...,    .2,    .1U286, 
.1,    .05. 
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FIG.  15>a.  See  caption  on  following  page . 
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FIG.  15a.     Graph  of  the  relationship  between  the  ratio    >f  -   (p-^Pq), 
where  p     is  the  pressure  behind  the  pure  gas  shocks  and  p     is  the 
pressure   in  front,   and  the   ratio   (u     0/ao)»  where  u     Q  is   (half)   the 
initial  discontinuity  in  the  transverse  velocity  and  a     is  the  sound 
speed  in  front  of  the   shock.     The  parameter  s^  is  the  square  ratio 
of  the  sound  speed  to  the   Alfven  disturbance   speed  in  front  of  the 
shock.     On  each  curve,  indexed  by  s   ,  the  point  corresponding  to  the 
minimum  value  of  u     n/a     at  which  cavitation  is  produced  is  indicated 
by  a  plus   (+)    sign.     The  dashed  curve  is  the   locus  of  these  points. 
The  point  on  each  of  the  incomplete  curves  corresponding  to  the 
"  critical"    values  of  u     Q/a     and  ^\ ,  -  that  is,  to  the  point  at 
which  the   switch-on  shock  in  the  resolution  reduces  to  a  pure  gas 
shock,  is  indicated  by  a  dot   (the  lowest  point  of  the  curves  indexed 
by  values  of  s     less  than  unity).     The  remaining  parts  of  these 
curves,  the  portions  below  the  critical  point,  are  shown  in  Figure 
15b.     The  scale  on  the  left  is  logarithmic. 
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FIG.  l$b.  See  caption  on  following  page. 
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FIG*  15b.     Graph  of  the   relationship  between  the  ratio  Vf  ■   (p  /p   )f 
where  p..    is  the  pressure  behind  the  switch-on  shocks  and  p     is  the 
pressure  in  front,  and  the  ratio  (u     n/a  )   of  (half)  the  initial 
discontinuity  in  the  transverse  velocity  to  the  sound  speed  a     in 
front  of  the   shocks.     The  parameter  s     is  the   square   ratio  of  the 
sound  speed  to  the  Alfve'n  disturbance  speed  in  front  of  the  switch- 
on  shocks^  the  presence  of  these  shocks  implies  that  s     is  less  than 
unity.     The   dot  at  the  upper  end  of  each  curve  corresponds  to  the 
"critical"    point  (u     0/a  ,  >i)   at  which  the  switch-on  shocks  reduce 
to  pure  gas  shocks. 

The  curves  depicted  here  should  be   regarded  as  the  con- 
tinuations of  the   incompleted  curves  of  Figure  15a.     The  upper 
portion  of  the    ,0625-curve  has  been  sketched  in  roughly.     The 
scale  on  the  left  is  logarithmic. 


See  caption  on  following  page, 
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FIG.  15c.  Graph  of  the  relationship  between  the  ratio  "*f  ■  (p-,/p0), 

where  p  is  the  pressure  behind  the  pure  gas  shocks  and  p  is  the 

pressure  in  front,  and  the  ratio  u  ~/b  .  where  u  _  is  (half)  th« 
r  '  y,0'  0*       y,0 

initial  discontinuity  in  the  transverse  velocity  and  b  is  the 
Alfve'n  disturbance  speed  in  front  of  the  shocks.  The  parameter 
8  is  the  square  ratio  of  the  sound  speed  to  the  Alfven  distur- 
bance speed  in  front  of  the  shocks.  These  curves  may  be  obtained 
from  the  corresponding  curves  of  Figure  15a  by  employing  the  rela- 
tion (u  _/b  )  ■  (u  _/a  )  \/s~  and  vice  versa.  The  scale  on  the 
y»0  o     y,0'  o'  v   o 

left  is  logarithmic o 
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FIG.  1U.  I:  Graph  of  the  relationship  between  vj    and  s  • 

Here,  vi    is  the  minimum  value  of  the  ratio  (p-,/p,J  of  the 
'   (cav  -i-  o 

pressure  p  behind  the  shock  to  the  pressure  p  in  front  of 
the  shock  at  which  cavitation  is  produced.  The  parameter  s 
is  the  square  ratio  of  the  sound  speed  to  the  Alfven  distur- 
bance speed  in  front  of  the  shock.  For  each  value  of  s  the 
corresponding  value  of  yi     is  obtained  from  the  logarithmic 
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FIG.  17a.  See  caption  on  following  page. 
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FIG.  17a.  Graph  of  the  relationship  between  the  ratio  (H  ?/H  ),  where 
—————  y  )5    x 

H  ,  is  the  transverse  magnetic  field  component  in  the  region  between 

the  rarefaction  waves  and  H  is  the  fixed,  longitudinal  magnetic  field 

component,  and  the  ratio   (u  0/a  ),  where  u  .  is  (half)  the  magnitude 

of  the  initial  discontinuity  in  the  transverse  velocity  and  a  is  the 

sound  speed  in  front  of  the  shocks.  The  parameter  s  is  the  square 

ratio  (a  /b  )  of  the  sound  a  to  the  Alfven  disturbance  speed  b^  in 
o  o  o  o 

front  of  the  shocks.  On  each  curve,  indexed  by  s  ,  the  point  corres- 
ponding to  the  minimum  value  of  u  n/a  at  which  cavitation  is  produced 

y  ,u     o 

is  indicated  by  a  plus   (+)    sign.     The  dashed  curve  is  the  locus  of  these 
points.     The  point  on  the  curves  indexed  by  values  of  s     less  than  unity, 
which  corresponds  to  the    "  critical*     values  of  u     0/aQ  and  ^    -  that  is, 
the  point  at  which  the   switch-on  shock  in  the  resolution  reduces  to  a 
pure  gas  shock,  is  indicated  by  a  dot.     It  is  easy  to  deduce  from  these 
curves  that,  for  fixed  values  of  u    Q  and  fixed  pressure  po     and  density 
p     in  front  of  the   shocks,  H     ,  is  a  non-decreasing  function  of  H» 
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FIG.  17b.  Graph  of  the  relationship  between  the  ratio  (H  -j/H„), 

—————  y>3  x 

where  H  ,  is  the  transverse  magnetic  field  component  in  the  region 
between  the  rarefaction  waves  and  H  is  the  fixed  longitudinal  mag- 
netic field  component,  and  the  ratio  (u  0/b  )  of  (half)  the  magni- 
tude of  the  initial  discontinuity  in  the  transverse  velocity  to  the 

Alf ve'n  disturbance  speed  b  in  front  of  the  shocks .  The  parameter, 

2 

s  is  the  square  ratio  (a  /b  )  of  the  sound  speed  a  to  the  Alfven 

disturbance  speed  b  in  the  undisturbed  medium  in  front  of  the  ad- 
vancing shocks.  These  curves  may  be  obtained  from  the  corresponding 
curves  of  Figure  17a  by  employing  the  relation  (u  nA>  )  ■ 

(u  _/a  )  >/s~  and  vice  versa.  It  is  easy  to  deduce  from  this 
v  y,0'  o  v  o 

figure  that  H  ,  varies  in  the  same  sense  as  the  pressure  p 
in  front  of  the  shocks  when  u  o5  the  longitudinal  magnetic 

y>o' 

field  H  ,  and  the  density  p  in  front  of  the  shocks  are  held 
x  o 

fixed . 
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Appendix  I 

The  Rantcine-Hugoniot  conditions  at  the  boundary  of  a  perfectly  conducting  neutral 
medium  which  expands  into  a  vacuum. 

Let   D(t)  be  an  arbitrary  time -dependent  domain  which  contains  part  of  the 
boundary,    S(t)   say,   that  separates  the  medium  from  free  space.     Let  E  (t)  be 
the  entire   surface  of  D(t).     Consider  the  following  system  of  equations: 

(a  )  f.n-0,  H"  •  n    -    0 

;r(t)  £j(t) 


(a2)  3F     I   ff  dV     -  H  U   «n  do"    -  n  x  E  d<T 

D(t)         £<t)  e  <*> 


(a  )  7  dV  +         J      ds     =      j  n  x  H  do~ 


D(t)  S(t) 


(I-D 


n(t) 


(a  )  w-  I       p  u  dV  -  p  u  (u«n  -  U«n)d<T 

>(t)  E  (t) 

-J         p  n  d<T    +  J  n   •  S(m)d6- 

E(t)  £<t) 


(a5)  It  J       p  dV     "    )  p(u-n-  U-n)  cKT. 

D(t)  E  (t) 

Here,  E"  and  H*  are   the  electric  and  magnetic  field  vectors,  7  the  volume  current 
density  within  the  medium,   J  the  surface  current  density  on  S(t),  p  the  density, 
p  the  pressure,  u  the  velocity  of  the  fluid,  n  the  outward  normal  to  E  (*)» 
tf'n  the  speed  of  the  surface  in  the  direction  n,  dC  the  differential  element  of 
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area  on  T~  (t)  and  ds  the   differential  element  of  area  on  S(t).     S         is  the 
magnetic   stress  tensor.     If  n  is  any  direction  in  space,  n«   S         is  given  by  the 
following  expression  * 

(1-2)  n.§(n,)  -  n(ff.n)  H  -  ^H2  n  , 

where  n  is  the  magnetic  specific  inductive   capacity  of  free  space.     Moreover, 
V  •   S         is  given  by  the  following  well-known  identity: 

(1-3)  7   •   "^   -  u.(v  x  H)  x  H  +  u-H  V   •   H   . 

In  the  medium  we  assume  that  the  electric  field  vector  is  given  by 

(I-U)  £"  -  n"H  x  u  . 

In  the  following  we  shall  suppose  that  dependent  variables  are  continuous 
and  have  continuous  space  and  time  derivatives  everywhere  except  possibly  across 
S(t).  Let  D  (t)  be  a  region  contained  entirely  within  the  medium.  Then,  the 
system  (i-l)  applies  in  D  (t)  if  we  set  J  ,  the  surface  current  density,  equal 
to  zero.   If  we  use  the  identities' 

It  f  Adv  =  j   §£  dv  ♦  J    a  u  .  H  d<r 


D(t)        r  (t) 


(1-5) 


n  *  Q  do"  , 

where  A  is  a  scalar  and  Q  a  scalar,  vector,  or  tensor  quantity  for  which  the 
operation  ■  #  "  makes  sense,  then,  on  taking  into  account  (1-3)  and  (i-u),  the 
system  (I-l)  reduces  to  the  following  set  of  partial  differential  equations  in 
D(t)l 


"5T 


"The  quantities  A  and  Q  are  assumed  to  be  continuous;  ^r  and  V  x  Q  are  assumed  to 
be  integrable  in  the  domain  over  which  the  integration  is  extended. 


-  112   - 


(b-j^) 


V   •  f  -  0, 


V   •   H  -  0 


(b2) 


3ff 
3T     + 


V  x   (H  x  u)   -  0 


(1-6) 


(b3) 


(V 


(b5) 


V  x  ff  -  J 


p(   |£  +  u    •   V  u)   +  Vp  +  (i(V  x  K)  x  H  ■  0 


3p 

at 


+  V    .    (pu)     -     0. 


Similarly,  if  D, (t)  is  a  region  contained  entirely  within  empty  space 
we  find: 


(c,) 


(1-7)       (c2) 


V  •  F  -  0, 


VxE+n|^     -     0 


V    •   H  -  0 


(c3) 


7  x  ff  -  0. 


Thus  the  system  (1-1)  is  equivalent  to  Lundquist  's  equations  within  the 
medium  and  to  the  quasistatic  approximation  to  Maxwell's  equations  in  free  space. 

Now  let  us  apply  (1-1)  to  D(t),  the  first  domain  discussed  in  this  Appendix. 
Let  S(t)  separate  D(t)  into  D  (t)  and  Tl(t),  D,(t)  being  the  region  in  free  space, 
Then  after  applying  (i-l)  to  D  (t)  and  D-,(t)  separately  and  subtracting  the  re- 
sulting equations  from  the  corresponding  equation  over  D(t),  we  get,  since  the 
volume  terms  cancel,  the  following  integral  relations  over  the  surface  S(t); 


(c^) 


ys(t) 


Hx[f] 


^     ds 


0,  n  x    |    |  ds     -     0 

'S(t) 


(1-8) 


(d2) 


(d3) 


'u  •  n   [H]   -  nx    [E  jlds     =     0 


'nx[H]-jlds     -     0 
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(1-8  cont'd) 


(d.  )  jp  u(u   •  n  -  U    •  n)jds  -   I       p  n  ds  -  n  .   S^ds 

'S(t)  "  S(t)  'S(t) 

(d  )  [p(u   .  n  -  U   •  n)     ds     -     0. 


S(t) 


Here,   the  normal  to   S(t),  n,   is  assumed  directed  into  D  (t);  and  if  Q  is  any 
quantity  (scalar  or  vector)    then 


(1-9)  [Q]     -     Q1  -Q0  . 

Since   D(t)  is  arbitrary,  £(t),  which  by  definition  is  that  part  of  the  boundary  of 
the  total  domain  contained  in  D(t),  may  be  chosen  arbitrarily  small.     From  (1-3) 
we  therefore  obtain  the  following  jump  relations  across  S(t): 

(e1)  [F]    •  n  -  0,  [  H  ]    •  n     -     0 

(e2)  [H]tr-H-nxCE]=0 

(ej  nx   [H]     =     J 

(1-10)         3 

(eu)  [mu]  -    [p]n  -  g  [  H2  ]  n  ♦  n(ff  •  n)    [H] 

(e   )  [>  3   x  °>     mi  "  PjC^Tj  •  n  -  U  •  n),         J     -    0,  1. 

Now  since  p1  vanishes  in  free   space  we  have,  using  (e^), 

(1-11)  n^     -     mQ     -     0 

It  follows  from  the  transverse  part  of   (e,  )  and  the   second  relation  of  (e-^  that 
£H  ]  vanishes*.     Consequently,   from  (e^ ,    (e2)   and  (ej  we  conclude  that    QE] 
and  J  vanish.     Moreover,  the  longitudinal  part  of  (e.  )  and    Qh  ]   -  0     implies 


It  is  assumed  here,  that  81    •  n  (■  H«    •  n)  does  not  vanish. 
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that    £'p  3   vanishes.     Hence,   since  p,-  0,  we   find  that  p     vanishes.     Since,  in 

addition,   p     is  related  to  p     by  means  of  the  adiabatic  equation  of  state  it 

follows  that  p     also  vanishes. 
o 

The  relations 

(1-12)  pQ  -  0   ,       P0  -  0  ,  QF]      -     0,  [IT]     -     0 

are  the  desired    "jump"    conditions  at  the  boundary  surface  S(t). 

It  should  be  observed  that  the  Rankine-Hugoniot  conditions  noted  in  the 
text  may  be  derived  in  a  similar  manner. 
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Appendix  II 

The  relation  of  q         to  q.      and  1-q         to  1-q.      as  q.      approaches  zero  and  unity 
cav  —    T,r  cav  T.r  —    T.r  -** — 

respectively. 

In  (U.13)   set  s  =  0  and  q  =  q       .We  then  have 

cav 


(II-D 


e(i-q     y 

cav 


(0+Dlog 


qtr(l*qcav)  1  1 


1qcav(1"qtr;  j  +  qcav  "  <*ti 


fa        k        1  U-V>9     "     <^%a,)9 


Let  a      and  hence  q         (  <  q.    )  approach  zero.     Neglecting  higher-order  terms, 
in  particular  the  logarithmic  and  sum  terms,  we  get 

(IM)  Jl.    .  q  <  JL.  -  -A.)  .  o. 

^tr  Hcav       ^tr 


Solving  for  q        we  have 
cav 


(n-3)  qcav  ~  &r  %r  ,   \ 


0. 


Now  let  q.   approach  unity.  Neglecting  the  first  two  terms  in  the  braces  in 
(1-1)  and  the  first  three  terms  in  the  summation  (it  is  assumed  that  ©  is  5)  we 
have,  after  suitably  rearranging  terms, 


(II-U) 


^cavY/2-5^-1" 


D*G£)-j-(*rc 


qt1 


0. 


Setting 

(II-5)               co 

1-q 
^cav 

we  then  have 
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(II-6) 


<o5  *   \  (l-q^  -   d-^) 


Assuming,   subject  to  later  verification,  that  ratio  of  the  first  to  the  second 
term  is  very  small  compared  to  unity,  i.e.,  that 
(3/2)(l-qtr) 


(H-7) 


GO 


<    <        1 


we  get 
(II -8) 


«>  -  (i-v) 


(3/2)^-1 


and  hence 


(II-9) 


GO 


<^> 


1/5 


(3/2)^-1 

\r 


1/5 


*ti 


1  . 


Clearly,  the  relation  (II-7)  is  satisfied.  Employing  (II-5),  we  have  finally 


(11-10) 


1  -  q    ~ 

cav 


d-V)U/5  (2/3)1/5-    V  -  1  . 
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